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ABSTRACT 


Analytical and experimental study for the effects of 
cooling system parameters on the heat transfer and 
temperature distribution in the electrode plates of a 
phosphoric acid fuel-cell has been conducted. 

An experimental set-up that simulates the operating 
conditions prevailing in a phosphoric-acid fuel-cell stack 
was designed and constructed. The set-up was then used to 
measure the overall heat transfer coefficient, the thermal 
contact resistance, and the electrode temperature 
distribution for two different cooling plate configurations. 

Two types of cooling plate configurations, serpentine 
and straight, were tested. Air, water, and oil were used 
as coolants. Measurements for the heat transfer coefficient 
and the thermal contact resistance were made for various 
flow rates ranging from 16 to 88 Kg/hr, and stack clamping 
pressure ranging from 0 to 3448 Kpa. 

The experimental results for the overall heat transfer 
coefficient were utilized to derive mathematical relations 
for the overall heat transfer coefficient as a function of 
stack clamping pressure and Reynolds number for the three 
coolants . 

The empirically derived formulas were incorporated in a 
previously developed computer program to predict electrodes 
temperature distribution and the performance of the stack 
cooling system. The results obtained were then compared 
with those available in the literature. The comparison 
showed maximum deviation of ± 11 percent. 
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CHAPTER I 
INTRODUCTION 


The fuel cell may be regarded as a continuously fed 
battery which converts the chemical energy of conventional 
fuel (such as natural gas and oxygen (from air) into 
electrical energy. Such a unit would not be reversible and 
cannot strictly be called a storage device; it is an energy 


con v 


ersion device or an electricity generator 


(Methane, I 
Methanol, I 
Naphtha etc.) 


Fuel 

Processor 


Fuel Cell 
jPower Section 


ower Conditioner 


Power 


Spent Fuel 


Figure 1 -Fuel Cell Power Plant Module 






Figure 1 is a block diagram of a phosphoric acid fuel 
cell (PAFC) power plant with three principal modules: the 
fuel processor, the fuel cell power section, and the power 
processor. 

Natural gas, methanol, and naphtha are principally 
considered for fuel cell use. Production of hydrogen, 
which is the major function of the fuel processor, occurs 
by reaction of the fuel with steam. 

The major components in the fuel processor subsystem 
are the reformer, two shift converters, and several heat 
exchangers. The reformer is basically a nonadiabatic , 
nonisothermal catalytic reactor which can operate as high 
as 1200 °C and 10 atm. — ■"" 

Within the fuel cell power section, hydrogen and 
oxygen react with a continuous production of DC electricity, 
waste heat, and steam (as a reaction product). The oxygen 
is obtained from air and the waste heat can be removed by 
the cooling system of the fuel cell stack. 

The heart of the power section is composed of a stack 
of several fuel-cell modules with four or five of such 
modules being sandwiched between two plates of the cooling 
system. The basic fuel-cell module can be made of two 
30 cm x 40 cm bipolar plates, which function as gas 
diffusion electrodes, and a phosphoric acid matrix which is 
sandwiched between them. 

The design and performance of the cooling system for 
the fuel-cell stack has significant effects on the 
operation and performance of the fuel-cell power modules. 
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In turn, th.se effects hs«e direct influence on the capital 
an „ operating cost of the fuel-cell power plant and the 

cost of electric energy produced. 

In addition to removal of heat generated by the 
olectrochemical reaction, the stack cooling system must he 
so designed that surfaces of the electrodes are kept at a 
reasonably uniform temperature of 200 "c for optimum 
operation and to prevent the production of thermal stresses 
that can cause distortion and subsequent failure of the 

electrode plates. 

The objective of this study is to analyse the heat 

transfer in the cooling system and to develop a correlation 

for the overall heat transfer coefficient in terms of the 

local leat transfer coef ficient . coolant thermophysical 

properties, fluid flow characteristics. thermal contact 

resistance and cooling channel configuration. 

!.! Qvgall Heat J r ansfer Coefficient In a Fuel Cell Stack _ 

Is some phosphoric acid fuel cell power plants, the 

stack i* arranged in such a way that each pair of cooling 

plates i sandwiched between five or more fuel cells modules. 

Si^le cell assemblies are stacked in a series 

bipolaTtrode. The directions of flow channels for air and 

hydro** fuel are perpendicular to each other . In addition 

to th,C€ power generated, heat is produced as a result of 

the comical reaction between hydrogen and oxygen. This 

heat jr removed from the cell stack by coolant such as (air, 

wateivetc...) passing through channels in cooling plates 

locakd approximately every fifth cell so as to maintain the 
locakd approxim OR^ALPJIpf 8 


p£roor«MLitt 
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'’igure 2 


is a schematic 


4 Cooling plate 


Fuel cell-fuel cell 
interface 


Cooling plate- fuel 
cell interface 


fuel cell stack may be 
le cell plate, convection 
:erface which is due to 
Lon through the cooling 
*een the channel wall and 
jut the overall heat 
all stack and its 
atical function will be 


discussed at a later chapter. 

1 • 2 Convective Heat Transfer Coefficienr 
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Engineering applications of convective heat transfer and 
mass transfer are extremely varied. In a multifluid 
exchanger we are concerned solely with the heat transfer 
rate between the fluids and the solid surfaces of the heat 
exchanger separating the fluids. Calculation of the 
temperature of a cooled turbine blade or the throat of a 
rocket nozzle involves convective heat transfer. The 
aerodynamic heating of a high-speed aircraft is a convective 
heat transfer process, but it also becomes a mass-transfer 
process when temperatures are so high that gas dissociates, 
forming mass concentration gradients. 

Obviously, theVcorabination heat transfer, mass transfer, 
and chemical reaction problem is the most challenging of the 
convective problems. Nevertheless, the bulk of this project 
is devoted to convective heat transfer. It is convenient in 
most cases to define a convective heat transfer conductance, 
or coefficient, such that the heat flux at the surface is 

the product of the conductance and a temperature difference. 
Thus : 

q " " h ’< T ! - T 2 > (1 ) 

The conductance h is essentially a fluid mechanic 
property of the system, whereas the temperature difference 
is, of course, a thermodynamic quantity. The usefulness of 

Eq ‘ (1) lies largely in the fact that in a great many 

* * * 
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technical applications it is close to being directly 
proportional to ( Tj - T 2 ), Tj and T 2 are temperatures of 
surfaces one and two respectively, as the linearity of the 
applicable differential equation reveals. Nevertheless, 
numerous nonlinear problems are encountered where h itself 
is a function of the temperature difference. It is important 
to note that this does not distroy the validity of Eq . (1) 

as a definition of h, although it may well reduce the 
usefulness of the conductance concept. 

^ ^ Ther mal Contac 

As aentioned previously in section 1.1, that contacts 
between the fuel-cell plate and the cooling plate in a fuel 
cell stack resulted in the existence of thermal contact 
resistance. This parameter is needed in the evaluation of 
the overall heat transfer coefficient for the system. 

Therefore, some historical background and definition about 
thermal contact resistance may be briefly discussed in the 
followiqj paragraphs. 

In the past, thermal resistances of contacts between 
solids reeived little attention, the contacts being assumed 
to be eiVer perfect or resulting in constant additional 
resistance. However, the quest for higher efficiencies 
required improved designs and materials pushed to the limits 
of their performance, and so the assumptions came to be 
increas&gl y questioned. This resulted in the investigation 
of the mriation of the thermal contact conductance with 
material properties, applied load, temperature, and fluid ) 
envird'eit. Hence, it is necessaTy^ to define thermal ^ 

ctHxmvmu# 
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contact resistance. 

1.3.1 Contact Resistance 

Whenever there is a flow of heat between two material, 
with thermal conductivities K a and Kj, , placed together 
under a load which is less than the load required for 
compressive yielding of the bulk material, actual contact 
between the surfaces occurs only at a discreet number of 
locations. This results in a temperature discontinuity 
across the interface. 

The lines of flow of heat converge in passing through a 
contact area which is smaller in size than the boundaries of 
the solids in contact. Thus, thermal contact resistance 
occurs because not all of the volume of the solids in 
contact is equally available for the conduction of flow of 
heat. There is an undirectional flow of heat away from the 
contact surface, but near the interface the flow becomes 
three-dimensional. The interface is a region of disturbance 
which extends into the specimen by approximately the distance 
separating two consecutive perfect contacts (54). 

From the basic Fourier equation, the linear temperature 
distribution some distance away from the interface is given 
by : 


£ 

k 


K a — 
a dx 



( 2 ) 


This may be extrapolated to the interface to determine 
the effective temperature drop ( AT ) across the interface, 

Figure 30 . If the thermal contact conductance h may be 

G0Mt&PJgpr$S 

pcmtjUAurtf 


defined by : 


h - aJA 

AT 


( 3 ) 


then 


h( AT) -Kil 
dx 


or h - K dt/dx 

AT 


The thermal contact resistance r is defined by: 


r c " 


1 

h 


( 4 ) 


Thus , if the heat flux and the temperature drop across the 
interface are known, one can determine the thermal contact 
resistance. It is important to realize that there is no 
localized interface contact resistance but rather a region 
of influence in the neighborhood of the contact. 


CHAPTER II 


REVIEW OF LITERATURE 

Estimation of the temperature profiles in an operating 
fuel cell is important for estimation of the power density 
distribution, thermal stability, and cooling requirements. 
Only a limited amount of information on this subject has 
been reported in the past. Baker and co-workers recognized 
this need and have performed a comprehensive steady state 
heat transfer in electrochemical systems (3,4,5). They 
studied various limiting and special cases to determine the 
maximum temperature of a stack; a two dimensional heat 
transfer analysis was carried out in the case of thick 
stacks where heat transfer in the direction of stacking was 
neglected. In the case of thin stacks, three dimensional 
heat transfer was considered with each wall at a different 
temperature. Infinite series solutions were developed for 
both thick and thin stacks.' The authors estimated the 
maximum stack temperature for the constant wall temperature 
case. An approximate formula to predict the effect of 
finite resistance to heat transfer at the wall, the effect 
of c«Ad or hot feeds, or non-uniformity of heat generation 
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was also carried out by using the method of Green's 
function . 

A single fuel cell with no lateral heat transfer and 
no conduction of heat through the cell in the direction 
perpendicular to the gas flow was considered (4). Heat 
transfer by conduction in the direction of the gas flow was 
considered negligible in comparison to the heat transfer by 
convection, and analytical expressions for the electrolyte, 
fuel and air temperature profiles were derived. 

Alkasab and Lu (6) have developed a heat transfer 
model for the phosphoric-acid fuel-cell stack combining 
mass, energy and electrochemical analysis. The derived 
three-dimensional mathematical model was utilized to 
develop a Fortran computer program in which this computer 
simulation included the determination of the effects on 
steady state temperature distribution in the cell plate and 
in the stacking directions: thermal conductivity, average 
current density, cell-plate dimensions, cell-plate size, 
coolant flow rate, cooling channels configurations, and 
inlet temperature of process air. 

Industries such as Westinghouse and Energy Research 
Corporation (45) have used a unique technique for cooling 
fuel cell stacks. One of the reactant gases, rather than a 
liquid, is used as the stack coolant. In the cooling gas 
scheme, the air feed stream is either split inside the 
manifold into both a reactant stream and a coolant stream 
which flow independently through the fuel cell stack and 
are either merged in the exit manifold or maintain separate 


reactant and cooling streams. 

Engelhard (46) has also developed a low cost cooling 
plate that can be made of conventional materials using 
conventional furnace brazing techniques. The coolant is a 
dielectric liquid. This cooler offers a potential cost 
reduction over other liquid and ebullient liquid coolers. 

A complete 5 Kw power plant incorporating a methanol 
reformer and a utility grid compatible power processor has 
been built and tested. 

United Technologies Corporation (47) has used an 
intercell cooling or two-phase water cooling method. Thin- 
walled, 2-pass copper tubes with stainless steel headers 
are the latest design. The acid environment requires the 
copper tubes to be coated with a thin teflon protective 
film; however, this somewhat reduces heat transfer. The 
chief disadvantage of these intercell coolers is their 
relatively high cost. The probability of cooler failure 
due to penetration of the film and subsequent corrosion of 
the copper has not yet been established. 

Table I (48) displays comparison of separate gas 
cooling to the two other commonly accepted cooling methods: 
process gas cooling and liquid cooling. 

The estimates are based on selected methods. No 
effort was made to optimize any of them; however, all 
systems are reasonable and are consistent for comparison 
purposes. 

Westinghouse has recently made a more detailed study 
of gas cooling. Under Doe Contract (49, 50) they developed 
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a lumped parameter fuel cell stack simulation code that 
calculates reactant gas compositon, current-voltage 
characteristics, and heat transfer characteristics for a 
g as-cooled fuel cell stack. In the model, the cell area is 
broken down to a grid of finite elements so that power and 
heat generation can be calculated as functions of 
temperature and reactant composition, as each varies from 

from point to point in a cell. 

Also, in the previous studies concern was primarily 

with the analytic evaluation of the conductance, or heat- 
transfer coefficient, under various conditions. The 
results of the evaluation were used, together with other 
heat transfer theory, for the analysis and design of fuel- 
cell systems. Some of these results have significant 
application outside the realm of fuel-cells, and there are 
cases where the conductance concept loses its useful 
significance. In these cases it will be more convenient to 
work directly with temperatures and heat transfer rates 

than to employ a conductance (41). 

Many authors (34) considered tubes with various flow 
cross section shape: a circular tube, rectangular tubes and 
triangular tubes. They considered heating and cooling from 
two surfaces and the effects of a peripheral heat flux 
variation around a tube. Next, they also considered a 
class of problems where the velocity profile is fully 
developed and remains fixed while the temperature profile 
is fully developed and remains fixed while the temperature 
profile develops. After this they became concerned with 


thermal-entry-length solutions for circular tubes and also, 
£ or rectangular tubes. A method was developed whereby the 
thermal-entry-length solutions for constant heat flux and 
constant surface temperature can be used to solve for the 
temperature distribution resulting from an arbitrary axial 
distribution of surface temperature or heat flux (34). 

The cooling channel configuration is important in 
analyzing and studying the overall performance of the fuel 
cell cooling system. However, until now emphasis has been 
primarily on analytic solutions, and experimental results 
have been referred to only where they serve to validate the 
assumptions used in building a mathematical model of the 
heat-transfer process. Nevertheless, experiments can and 
do form a primary source of convection heat-transfer data 
for engineering applications; if the flow geometry is 
complex, it is often far easier to perform experiments than 
to attempt to deduce heat-transfer rates by analysis. Since 
the present work deals mainly with the convective heat 
transfer analysis, it becomes essential to investigate and 
acquire an experimental technique to measure the heat 
transfer coefficient where knowledge of it in engineering 
is required for a very wide range of activities. Numerous 
straight-forward and ingenious techniques have been used to 
determine the convective heat transfer coefficient 
experimentally. Moreover, in general, the scope of this 
review is confined to techniques for determining the steady 
state heat transfer coefficient and to summarize the 
available techniques along with a selection of their 


Hoar Tran sfer Coefficient. 

Baughn (80 has recently summarised some methods of 
gisectly measuring hsaf transfer for forced convection fa 
Butts, so she examples chosen for each -hnique are for 
flat surfaces. Davenport (9) used a straight forvar 
technique for the mean heat transfer coefficient for the 
air-side of loeered surfaces in eater cooled heat exchanger 

specimens. It was simply based on evaluation of the 

j mas c flow and temperature 
enthalpy change from the measured mass 

. rhe test surface, 

change of a liquid on the other si 

Ma „ y authors considered electrical heating, which 
rhe most common experimental technique of determining loca 

and mean heat transfer coefficients. The electrical power 

Ki»i- rhe energy losses have 
• is easy to measure but the energy 

energy input is easy 

to be carefully evaluated. 

Three techniques have recently been used. Kim et 

al. (10) used 12 copper plates 460 mm wide. 50 mm long an 
6 mm thick for a stud, of full coverage film cooling in a 
wind tunnel. The plates were heated b, resistance wires 
installed in slots machined into the back side of eac 
p late . To minimise the energy loss from the back of the 

plates, heated water tubes were used. 

a a a nil used 18 stainless steel 
Eriksen and Goldstein (ID usea 

„ wid . 50 mm long and 0.25 mm thick for 
foil heaters 203 mm wide, 

' dies in a wind tunnel. The foil was 
film cooling studies in a 

v uesrf with 50 mm of styrofoam, and 
electrically heated and backed with 
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had thermocouples on its surface. Blair (12), 

Hand, obtained a acre detailed knowledge of tbe beat 
transfer coefficient distribution b, -sing U* separate 
, «, thick composite heaters on an area 508 «« * 122° "" 
again for film cooling studies in a wind tunnel. The 
heaters were backed with about 15 mm of rigid urethane foam, 
a„d the infrared technics was used to measure the wa 

temperature (13). 

Another technique is to use electrically conducing 
vail coatings. Two types are commercially availble. a 
transparent type of vapor -deposited gold on polyester film 
0 2 mm thich * and a carbon impregnated coating on a 

plastic sheet 0.31 mm thick. These hare been successfully 

** frtr a variety of flat 
nsed with liquid crystal sheets for a 

and curved geometries by Hippensteele (H). This 

combination gives local qualitative and quantitative data. 

Heat flux sensors -ere used with the Fourier heat- 
cnnduction equation to obtain local heat transfer 
coefficients. It is important that the, do not disrupt t e 

11 ,. alter its temperature.. Crawford 

gas flow at the wall or alter its t p 

6 n flat film-cooled 

at al. (23) have used foil-type sensors on 


♦From Sierracin/ Index 

**Such as those from 


Products, Chatworth^ CA, 
Liquid Crystal Technology 


USA. 

LTD, 


LONDON, UK. 
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surfaces, while Hay and West (24) have constructed and 
applied slug-type sensors to swirling flow in a pipe. 
Miniature sensors which are about 25 times smaller, x.e. 

0.7 mm x 0.2 mm x 0.8 micrometer thick, have been developed 
by Portrat, et al. (17) for turbomachinery. 

Other methods of indirectly measuring heat transfer 
were conducted. For example, mass transfer experiments are 
easier to set-up, have cleaner boundary conditions, are 
easier to study, and are more accurate than the 
corresponding heat transfer process. The heat and mass 
transfer analogy approach where the gas concentration at 
the walls plays the same role in the mass transfer process 
as does the wall temperature in the heat transfer process, 

i . e . : 

h (T u - T„ ) - h m <C u - C_ ) (5) 

The simplest and most useful form of the analogy was 
proposed by Chilton and Colburn (55). 

This is : 

St . (Pf) 2/3 - St m . <Sc) 2/3 (6) 

and where the flow fields are identical for. the heat and 
mass transfer systems! 

2/3 


h - h C (Sc/Pr) 
m p 


(7) 



So, by measuring h ffl , h can be determined. 

A comprehensive account of the analogies between mass, 
heat and momentum transfer is given by Sherwood et al. (56). 

2 . 2 Thermal Contact Resis tance 

Jacob and Starr (21) investigated thermal contact 
resistance of interface joints between various metals in 
vacum as a function of pressure at room temperature and at 

the temperature of boiling nitrogen 93 K. 

Brunot and Buckland (22) investigated the dependence 
of thermal contact resistance on laminated and cold rolled 
steel joints of various surface roughness and -contact 
pressures. As expected, contact resistance decreases with 
increased pressure. Their test on solid steel blocks with 
various degrees of smoothness indicate that the smoother 
the surface, the lower the resistance. They infer that it 
is probably " due to the thinner layer of air or the larger 
area of contact encountered ". Another conclusion is that 
tests on laminated steel blocks indicate that if a thick 
metal shim is used between the two surfaces, the hardness 
of the shim has little effect. If aluminum foil is 
substituted, the resistance is lowered as pressure is 
increased. They ascribe this to the fact that the 
laminations are embedded deeper in the foil, so that there 
is more effective conformity between surfaces. The effect 
of interstitial fluid on thermal contact resistance was not 

measured . 

Weillis and Ryder (23) determined the dependence of 
thermal contact resistance of joints on pressure, surface 


finish, temperature, interstitial fluid, heat flow, and 
temperature drop. Their results show a linear dependence 
of thermal contact conductance of dry steel joints and an 
exponential dependence for aluminum and bronze. Other 
significant results are the decrease in thermal resistance 
with a decrease in surface roughness of both dry and oil- 
filled joints, and a higher resistance of dry joints to 
oil-filled joints. The effect of the oil decreases at 
higher pressures. The thermal resistance -is decreased by 
c o ppe r— p 1 a t i n g on the surface of a joint. 

Barzelay et al. (24) investigated the effect of heat 
flow, the temperature of the joint and later also 
investigated (25) the effect of pressure and dissimilar 
metal combinations on thermal contact resistance. They 
found that for the same metal combination the thermal 
contact resistance depended upon the direction of heat 
flow, whether heat flowed from steel to aluminum or from 
aluminum to steel. They found that, other factors 
remaining the same, the thermal contact conductance for 
aluminum to aluminum contact was greater than that for 
aluminum to steel contact, which in turn was greater than 
for steel to steel contact. The least value was observed 
for steel to aluminum contact. The reason for this was not 
satisfactorily explained by the authors, as pointed by 
Powell et al. (26). 

They also found that interface conductance increased 
with pressure and with the mean temperature at the interface. 
Another conclusion was that surface roughness alone was not 


a dominant parameter in determining thermal conductance of 
contacts; overall flatness had a more important role in 
determining the configuration of surface matching. 

It appears that the majority of researchers have not 
included the effects of this parameter (contact resistance) 
in their analysis that involved fuel cells. Also, in the 
evaluation of the overall heat transfer coefficient, the 
terms " contact pressure " and " contact resistance " were 
not considered either. More details about the approach 
used in measuring the contact resistance and the 
significance of it will take place in a later chapter. 


CHAPTER III 


SYSTEM DESCRIPTION 


. In this chapter, a general description of the apparatus 
and the test procedure are presented. Also, shematic views 
for devices used are illustrated, and measuring instruments 
and other equipment are described. 

3 . 1 Construction of the Fuel Cell 


The objective of this task is to design and construct 
an experimental set-up that simulates the output and 
working conditions that prevail in the fuel-cell module and 
the cooling system of a phosphoric acid fuel-cell stack. 

The set-up is assumed to simulate a fuel-cell module 
operating at an average temperature of about 190 °C and one j 

i 

atmospheric pressure, having an ouput of 0.8 volt and an 
average current density of about 0.325 amp/cm’. 

The fuel-cell module shown in Figure 3 consists of 
two 0.30 m x 0.41 m mating graphite plates with an upper 
plate thickness of 6.35 mm and a lower one of 9.5 mm. 

Nichrome wire was used to hold the heating loops as shown 
in Figure 4 

As was stated before, graphite is the principal 
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CELL PLATE 
(upper purr) 


MICA SHEET 

1.59 nun HOLE 


MARIN (TE 

HEATING ELEMENTS 



FIGURE 3 - ILLUSTRATIVE CELL CONFIGURATIONS 


1. MARIN fTE I 

2. NICROME WIRE 

3. ELECTRICAL CONNECTIONS 


Figure 4 


- Heating Elements Configuration 



material used to construct fuel cells. Contact with 
nichrome heating elements caused circuit shortage because 
the graphite is an electrically conductive material; 
therefore, the use of marinite became necessary in order to 
electrically insulate the heating elements. 

Marinite is structural insulation not affected by 
moisture or high humidity, and will not rust or corrode. 

In addition, it withstands a temperature of up to 650 °C . 

A sheet of marinite 0.25 x O.SSm’with a 6.35 mm thickness was 
sandwiched between the two plates, after milling 4 circular 
loops of grooves 4.76 mm in diameter in each plate to fit 
the heating elements. This sheet of marinite was inserted 
into the lower plate of the cell module covering an area 
of 150 inches square. The face with circular grooves which 
contained the heating elements was in touch with the bottom 
part of the lower plate. In order to prevent electrical 
contact between the nichrome wire and the upper plate, a 
very thin sheet of mica was inserted in between. Thus, the 
heating elements, the marinite and the mica were sandwiched 
between the two graphite fuel-cell plates and bolted 
together, forming the simulated fuel cell module. 

3 . 2 The Experimental Set-up 

A general view of the test installation is shown in 
Figure 5 , a top view is shematically illustrated in 

Figures 6,7,8, and a close-up of the pressure 
loading system showing the location of the test assembly 
is shown in the shematic of Figure 9 • 

The set-up consists of the following: 



Figure 15 - The Experimental Set-up (Air Case Shown) 



1 TANK FOR COOLANT (OIL) 

2 RUMP 

3 FLOW METER 

4 FUEL CELL PLATE 

• COOUNG PLATE (SERPENTINE TYPE 
SHOWN) 


• THERMOCOUPLES CONNECTED TO DATA LOGGER 
7 FLUKE OATA LOGGER 

• OIL COIL ICE BATH 

• WIRE CONNECTIONS TO HEATING ELEMENTS 
10 POWER SUPPLY 


Figure 6 * Flow Diagram of the Apparatus when Oil is Used as a Coolant 
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1 BLOWER 

2 PITOT TUBE 

3 MANUFOLD 

4 FUEL CELL PLATE 


3 COOLING PLATE (STRAIGHT CONFIGURATION SHOWN) 

6 THERMOCOUPLES CONNECTED TO DATA LOGGER 

7 FLUKE DATA LOGGER 

8 WIRE CONNECTIONS TO HEATING ELEMENTS 

9 POWER SUPPLY 


Figure 7 - Flow Diagram of the Apparatus when Air is Used as a Coolant 





1. TANK FOR COOLANT (WATER) 

2. PUMP 

3 . FLOW METER 

4. CMC UNIVERSAL PRESET COUNTER 

5. FLUKE DATA LOGGER 
4. FUEL CELL 


7. COOUNG PLATE (SERPENTINE TYPE IS SHOWN) 

8. THERMOCOUPLES CONNECTED TO DATA LOGGER 

9. DRAINER 

10. WIRE CONNECTIONS TO HEATING ELEMENTS 
tl. POWER SUPPLY 


Figure 8 


- Flow Diagram of the Apparatus when Water is Used as a Coolant 








1- Fuel-cell module 

2- Fuel-cell cooling plates 

3- Power supply unit 

4- Temperature measurement and recording equipment and 
instruments . 

5- Coolant circulations and fluid flow measuring 
devices . 

6- Stack pressure loading system 

3.2.1 Fuel-Cell Module 

Heat generation in the fuel— cell module was implemented 

\ 

by supplying the required amount of wattage to simulate an 
actual operating fuel cell from the power supply. The heat 
source consisted of 4 loops of coils of 80-20 nickel- 
chromium alloy resistance wire. These coils were centrally 
positioned inside machined circular grooves in the marinite 
sheet with the ends of each coil loop connected to the 
power supply. The terminals of the coil were electrically 
insulated by compacted magnesium oxide which has the unique 
quality of being a good insulator at high temperature, and 
at the same time a good heat conductor Figure 4 . Each 

heating loop was connected to an autovariac which supplies 
it with 120 volt A . C . By varying the voltage across each 
loop, the shape of isotherms across any section in the heat 
assembly also could be varied. This gave an additional 
control in order to get a uniform heat flux across the cross 
sectional area of the fuel cell plate. 

3.2.2 Cooling Plates 

The cooling system consisted of two types of cooling 
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plate configurations and three different types of coolants. 
The serpentine configuration was used for oil and water and 
the straight channel configuration for air. 

A . Serpentine Configuration 

As shown in Figure 14 , this configuration consisted 
of two mating plates of graphite, size 0.30 x 0.41 m 1 , and 
a copper tube 9.5 mm in diameter, 3.66 m long. The copper 
tubing was inserted as loops inside the cooling plate. 
Grooves were machined in each half of the cooling plate to 
fit four and a half loops of the copper tubing. The 
coolant was circulated through the system by the use of a 
pump with a variable speed motor and a maximum capacity of 
5 gpm. The variable speed motor provided control of the 
coo'laat flow and the ability for it to operate at different 
flow rates. 

3 . Straight Channel Configuration 

The straight channel configuration is shown in 
Figures 15 , 16 . This cooling plate was built from 

graphite of the same dimensions as that of the serpentine 
type- Two mating halves were constructed with each half 
haviig a total of 30 half channels; the combination of the 
two ialves formed the straight channel configuration. 

In this configuration, air was forced through the 
manifdlL system by a blower, a 51 mm diameter plastic tube 
was u^d' to connect the blower outlet to the manifold 
system inlet. To minimize corrosion, the manifold was 
fabricated from steel and the tube was constructed from 
h«a vwr duty plastic. 


OffidNALPMt tt 
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3.2.3 Power Control Unit 

This unit consisted of the metal panel where a set of 
6 variacs, ammeters and voltmeters are dispersed in series 
of 3. Each independent set of one (variac, voltmeter, 
ammeter) was connected to the heating loop by a heavy duty 
electric wire. Monitoring the amount of electric heat in 
each loop was easily managed and the amount of power was 
measured with considerable accuracy. This instrument 
obtained its electric power from the main switch that is 
provided by the university in the lab. 

3.2.4 Temperature Recordings Devices and Instrumentations 

In this experiment, measurement of temperature was the 

most important factor and the most sensitive aspect as far 
as accuracy is concerned. Thermocouples type T copper- 
constantan were used for measuring the temperature. Over 
100 temperature measurements were needed during the 
operation. Therefore, a data acquisition system was 
installed, type Fluke (Model 2200 B), along with an 
extender chassis (Model 2201 A) giving the system a 
capacity of accomodating 100 channels. Basically, this 
unit is designed to scan and measure analog input data, and 
then output the measurement data in digital form. Output 
data is displayed and can be recorded. The control section 
operates under program control to select input channels, 
and scale the resultant measurement readings for display of 
recording purposes. 

A total of 72 thermocouples Figure 13 , were installed 

to measure the surface temperature of the cell plate. 

7 ^ * r * 
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Four thermocouples, which measured the wall temperature, 
were attached on the plate centerline along the length and 
at four locations 0.15m apart axially from each other. The 
temperature for inlet and outlet of coolant were also 


} 


measured and monitored by the Data j^ogger ♦ Thermal contact 
resistance temperatures were recorded. Other temperature 
measurements such as room and surroundings were obtained by 
a Fluke Digital thermometer. The thermocouples were a 24 
Gauge or0.25mm in diameter, which were glued to the lower 
cell plate through holes 1.59mm in diameter made on the 
upper plate, and then were connected to the Data Logger. 
Also, the rest of thermocouples were glued to their 
locations using Viton cement which not only holds the 
thermocouples in place but withstands high temperature and 
provides good heat conduction. 

3. 2 . 5 Coolant Circulation and Flow Rate Equipment and 
Instruments 

Flow rate for the cooling air was determined by using 
a pitot tube. The flow is manually controlled by varying 
the speed of the blower motor; this arrangement was done by 
using a variac to draw different voltages across the blower 
motor in order to obtain the desired flow rate. 

Water flow rate was measured by using a turbine flow ^ 


transducer, an inline m elf eTiifg~~de vice utilizing a bladed 
rotor to generate flow information. It has the ability to' 
measure liquid and gas under high temperature and pressure 
conditions with accuracy and reliability. It was connected 
to a flow counter CMC (Model 614 A) where a digital readout 
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Figure 12 — View for the Fuel Cell Module 




of the frequency is obtained, and a calibration curve of 
f rfitjusncy versus gallon/minute was established, A 
mechanical flow meter (Model 1110 B) Brooks with a maximum 
capacity - of 0.40 gpm was usedf to determ ine the flow rate of 
oil. A certificate of calibration was provided^ by_ the 
manufacturer and a manual calibration was imposed to check 
over the data. Water and oil flow rates were manually 
controlled by a special valve attached to the pump used in 
this experiment, 

3. 2 . 6 Stack Loading System 

This system simulated the pressure applied to the fuel 
cell stack to hold the fuel— cell module and the cooling 
plate together. As shown in Figure 12 , the stack of the 

fuel-cell module and the cooling plate were sandwiched 
between two steel plates having a thickness of 12.71 mm each. 
Pressure was applied to the stack by a 4 ton capacity 
hydraulic jack and measured by gauge as shown in Figure 9 . 
3. 3 Test Procedure 

A shematic of the flow system is shown in Figure 6 
The test began by applying electrical power to the heating 
elements in the fuel cell module. Temperatures at 
different locations in the cell plate were monitored 
until a uniform heat generation and an average cell 
temperature of 190 °C were obtained. At this point, the 

cooling system was put into operation with coolant using 
oil, water and air as coolants. The oil and water were 
circulated through the serpentine passages or the air was 
forced through the cooling channels. 
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When oil was used as a coolant, a cooling coil, which 

is shown in Figure 6 • has been designed and built so that 

* 

the oil is cooled to 21 °C which is the oil inlet 
temperature used in this experiment. Temperatures were 
monitored again until steady state conditions were achieved. 

Tests were conducted for various interface pressure (0 to 
3448 Kpa ) . 

Temperatures were digitally recorded and were 
continuously available on a 5 minute update cycle. The S / 

control of temperature of the coolant was achieved by an 

1 / 

isothermal bath (oil case); the copper coil was contained in 
a 50 liter capacity reservoir tank. The bath temperature 
was controlled with simultanuous hot and cold water. The 
temperature was monitored by the Data Logger during the test 
period providing continuous reading for the inlet 
temperature. Bath temperature was sensed by a thermocouple 
installed in conjunction with the experimental thermocouple 
network. Water temperature was pumped from the reservoir 
t a n k while the temperature monitored by the Data Logger and 
controlled by a similar method as oil achieved the desired 
temperature. Air temperature was maintained between 20 °C 
and 23 °C during all the test period, an average of 21 °C 
was used in carrying out the mathematical calculations. The 
wall temperatures were measured at 4 stations 15.25 cm axially 
apart from each other. An arithmetical average for all four 
was used in the calculation. Twelve thermocouples were 
installed to measure the temperature drop across the 
interface: three were stationed at 3 locations 6.33 cm axially 


apart from each other, 2 sets in the cell plate and the 

I 

other 2 in the cooling plate. The effective temperature 
drop was found by linear extrapolation to the interface as 
shown in Figure 30 

The pressure was measured by a zero to 6895 Kpa gauge; 
the gauge had been calibrated by the general instrument 
corporation with a piston gauge which has been compared with 
master piston gauges. Furthermore, the calibration 
certificate suggests an accuracy of 0.1 percent for the 
gauge employed in the program ; 

The heat was supplied by the loops of nichrome wire, <" 
gauge 30 of about 8.8 ohms at room temperature. The 
electrical input power was determined from the current and 
voltage measurements. All measuring instruments in this 
part were calibrated and carried certifications of the test 
showing a maximum deviation from nominal of 0.01 percent. 

The coolant flow rates were measured by flow meters 
except for air where a pitot tube was employed , power 
input was constant during the test period. Conditions 
were then changed by varying one of a number of parameters 
which included the coolant flow rate, the interface pressure, 
as well as the cooling plate configurations. Each test was 
repeated three or four times and measurements were recorded 
accordingly . 

For completness, the experimental work included careful 
insulation of the (cell/cooling) plate to minimize heat 
losses. Non-asbestos millboard type Binder-Cement that 
withstands temperatures of up to 538 °C was used for this 


purpose. In-order to assess the magnitude of the heat loss, 
the apparatus .as operated' it lo. power .with no coolant flow 
and the wall temperatu re wai mias ured. These data were then 
used to jTorrect The 

total heat loss was^found to be approximately 1 X of the 
heat input. The corrected power input was employed in the 
evaluation of the average fluid teip'erlture-'gradient”. Data 
reduction followed the normal pattern of the evaluation of 
Reynolds numbers, thermal contact^ resistance, and finally 
Sverige heat t ransf er_co .f f ic ten t 7I7T.~ hJIT transfer 
coefficient based on average surf ace Ti'TTTfliVTnd average 
wall temperature to the average fluid temperature 
difference) . 
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FIGURE 16 - SECTION A- A AS 8HOWN FOR STRAIGHT CHANNELS CONFIGURATION 
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CHAPTER IV 
ACCURACY OF DATA 

Since there could be many sources of error with 
varying degrees of influence for different sets of readings, 
it is therefore appropriate to discuss individually these 
sources of error and their influence on the final result. 

4 . 1 Errors Due to Thermocouple Readings 

Any error in the thermocouple readings is of great 
significance as it is from this, directly or indirectly, 
that all the numerical results were obtained. Each of the 
thermocouples used in these tests were made of wire from 
the same lot and read on the same Data Logger. 

Obviously, the significance of this error would depend, 
on the temperature drop across the interface which is used 
in the evaluation of the thermal contact resistance. 

However, the accuracy of measuring the overall heat 
transfer coefficient depends on the accuracy of the thermal 
contact resistance values. At high pressure the minimum of 
temperature drop across the interface is found to be 2 °C 
as seen in Figure 30 . When scales were used, the 
temperature drop across the interface varied from 2 °C to 
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6 °C . Inlet temperature of the coolant was between 20 °C 
and 22 °C for all types of coolants. Obviously, the 
maximum error, due to thermocouple readings, in presence of 
the Data Logger was not more than 1 percent. 

4 . 2 Errors Due to Heat Leaks 

There are possible paths through which the surroundings 
can exchange heat energy with the test module where the 
essential measurements are taking place. 

The first is by convection through the body of the 
test module. To eliminate the convective component of heat 
transfer to the surroundings, the fuel cell module has been 
insulated. The quality and the thickness of the insulation 
used were such that heat loss to the surroundings was 
expected to be insignificant. It was calculated to be 
about 15 watts which is less than 2 percent of the total 
heat generation. This could result in some percentage of 
error in the calculation of the heat transfer coefficient 
and conductance values. 

The second source of error was in the uncertainty of 
the uniformity of the heat flow which could be caused by: 

a. Heat loss 

b. Non-uniformity in the heat path 

The non-uniformity of heat flow has two origins: the 
first is the non-uniformity inherent in the very nature of 
contact resistance; the second is the disturbance created 
by the thermocouple insertion. Nothing can be done with 
the first one. The second one was minimized by spot- 
welding the tip of the thermocouple to a small brass 



sphere 1.27 mm in* diameter. This is done so that an average 
temperature over a small region at the center could be 
measured. In this way any inaccuracies due to thermocouple 
bead not being exactly at the center were diminished. 
Moreover, in order to have uniform horizontal isotherms in 
the cell plate itself, the latter must have boundary 
conditions. This was achieved by the use of graphite in 
the construction of both the cell plate and the cooling 
plate. Besides, this material was chosen because its 
electrical resistivity did not vary with temperature, 
resulting in a uniform heat flux. 

Precise measurements of the axial location of 
thermocouple was necessary for the accurate calculation of 
the temperature gradient. This is especially important in 
the determination of the temperature drop across the 
interface. 

Although thermocouples were placed as close as possible 
to the interface, a linear extrapolation of the axial 
temperature gradient could have been inaccurate when the 
temperature drop was unusually small or when a radial 
gradient of sufficient magnitude was present. However, any 
error due to this extrapolation was assumed negligible. 

^ • 3 Errors As sociated with Observation of the Properties 

to be measured 

This section is a brief discussion of errors made due 
to errors in calibration and readability of pressure gauge, 
flow meters, voltmeter and ammeters. 

As indicated earlier, the pressure gauge has an 
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accuracy of 0.1 percent. This means that for a range of 
maximum pressure of 6895 Kpa used in the present research, 
the pressure gauge had an error ±6.9 Kpa. To account for 
errors in readability and calibration we shall adopt a 
maximum of ±35 Kpa. 

The test flow meters have an accuracy of 0.05 percent 
for the turbine flow transducer and 0.1 percent for the 
Brooks instrument type. The error due to measurements of 
flow from the pitot tube is accounted to be about 1 percent 
Thus, the maximum combined error due to flow rates 
measurements would amount to 0.4 percent. Other 
measurements such as readings of voltage, current and heat 
losses along thermocouple wires, were assumed to be about 
one percent. 

Nevertheless, the error in calculating the heat 
transfer coefficient for 3 different tests under the same 
conditions for oil, water and air was found to be less than 
one percent. However, all the above factors might have 
influenced the values presented for heat transfer 
coefficient and thermal contact resistance to an estimated 
± 10 percent in the extreme. 
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The energy balance' for the coolant includes convective 
heat transfer between coolant and channel walls and con- 
vection by the coolant. The energy balance for the coolant, 
for a unit length along the channel, can be written as 


hs <T W - T c ) - .C o dT. . 0 


dx 


where 


h ■ Heat Transfer Coefficient 

w/m* . °C 

S ■ Perimeter of Cooling Channel 

m 

Ig* Average Coolant Temperature 

°C 

Average Wall Temperature 

°C 

m ■ Cooling Mass Flow Rate 

Kg/hr 

C p - Specific Heat of the Coolant 

KJ/Kg . °C 

x » Distance from the Edge 

m 


The boundary condition is: 

T c * , at inlet at x « 0 


In carrying out this mathematical analysis the 
following assumptions were made: 

(a) Thermal conductivities, specific heats, and densities 

® ^ ^ fluids involved are assumed to be constant. 

(b) The heat flux at the wall is constant. 

(c) The edges of fuel cell plate and cooling plate are 
assumed to be adiabatic. 

(d) The coolant flow is assumed to be steady. 

Thus, if the flow enters at temperature T^ , the 
solution from Eq . (8) yields a value of the exit 
" temperature, T q , in the following form. 
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To " T i + < T w - T t )( 1 - e -<5) 


(9) 


where , 


d(x) . hSx 
o . C 


( 10 ) 


Eq. (9) was solved numerically on the IBM 370 computer 
for five different flow rates and for air, oil and water 
properties. Comparison of both the calculated and the 
measured values for the exit coolant temperature are 
presented in Appendix (B). 

5,2 Heat Transfer Analysis for the Ser p entine Configuration 

In tuba flow, the heat flux q" is usually defined by 
the relation. 


or , 


and , 


" c p c T b 2 - V 


V ■ h ( T w - T b ) 


* <T W - T„) 


( 11 ) 

02 ) 

(13) 


where 


m 

°C 


A - local Surface Area for Heat Transfer 
T w * Garage Wall Temperature 
T b“ * Bra 8e Fluid or Bulk Temperature «C 

Subsctpsts 1 and 2 refer to inlet end outlet conditions 

In leniving the mathematical relation for the convective 


<*pom vSu* 
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heat transfer coefficient in the serpentine cooling plate 
configuration, the followings were led to be : 

1- Wall temperature is the logarithmic average of the 
measured temperatures at various locations at the wall. 

2- T b , is the logarithmic average of the inlet and 
oulet temperatures of the coolant. 

The heat transfer rate through the surface of the 
cooling plate was determined from the measurements of the 
electrical power dissipated by the nichcrome heating 
elements in the cell plate, the outer wall temperature, and 
the mean of the cooling fluid temperature as expressed in 
Eq. (13?. 

The electrical power was determined from the heating 
elements current and voltage measurements. The inaccuracy 
in the current and voltage measurements was estimated to be 
less chan 5 percent. This is believed to be the major 
uncertainty in' the heat transfer measurements. Check 
me a s tneaen t s were made of the voltage drop across the 
heating loops . This gave an upper limit of the power 
generated in the cell plate, since it included losses due 
to coirtect resistance. Measurements showed a little 
variat i>jr in the electrical resistance per heating loop 
over tJe surface of the cell plate. 

Tie outer wall temperature was measured with 
thermocouples at various locations as shown in Figure 19 
Inlet tjemperature of coolant was measured with a 
themocouple located in a well mixed location at the 
entnace to the test section. The outlet temperature was 
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3 



2 - Fuel-Cell Plate 

3 - Thermocouple Location 

At The Wall 


Figure 18 -Measurements of Wall Temperature 

measured in the settling chamber at the manifold in the air 
case and at a probed thermocouple inside the copper tube in 
the oil and water cases. 

Reynolds and Nusselt numbers were determined from the 
following relations. 



where , 

V ■ Velocity 

D • Cooling Channel Diameter 


m/ sec 
m 
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vr « Dynamic Viscosity of the Fluid m’/sec 


NU - h D 

K 


( 15 ) 


where 

K » Thermal Conductivity of Coolant w/m.°C 

♦ 

For oil and water Eq . (13) was used in the calculations of 
the heat transfer coefficient. 

5 . 3 Heat Transfer Analysis for the Strai ght Channel 
Conf iauration 

The heat transfer coefficient for the cooling air was 
determined from the energy balance between the heat 
removed by the coolant and the convective heat transfer 
betweei the channel walls and the coolant. 


Figure ft- - Geometry of Cooling Channel (Straight) 

At steady state the energy balance can be written as 
follows: 


M (T w - T c ^ 


“• C p < T o 


T i > 


(16) 


and 


"■• C n( T „ - T 1 l 

* <T W - T c ) 


In. 


(17) 
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where 


m - Cooling Air Flow Rate Kg/hr 

A - Heat Transfer Area for the Cooling Channel. m 

T * Outlet Temperature of Air 

0 

T » Inlet Temperature of Air 

1 

T - Average of Inlet and Outlet Temperature °C 

of Air 

T = Wall Temperature 


°C 

°C 


If 


¥ 


In this analysis the hydraulic diameter is used for 
evaluating Reynolds number. This is so because most of the 
velocity change in a turbulent flow takes place very close 
to the wall surface and is relatively independent of the 
proximity of other wall surfaces. For this reason the 
shape of the flow tube cross section has little effect on 
the shear stress at the wall, except where sharp corners 
are involved. The hydraulic diameter is defined by: 


'H 


2 x y 


( 18 ) 


x + y 


x - Width of Cooling Channel 
y - Depth of Cooling Channel 


m 

m 


.5.4 The Overall Heat Transfer Coef ficient 

The heat-transfer process may be represented by the 
network shown in Figure 20 , and the overall heat transfer 
is calculated as the ratio of the overall temperature 
difference to the sum of the thermal resistances. 





r 

c 



KA 


1 _ 

hA 


Figure 20 -Overall Heat Transfer through the ( Cel 1 /Cool ing ) 
Plate Unit. 

This is represented by the following relation, 

T s - T c 

q - — — • (19) 

1 + Axj + ax 2 + r c 

hA KA“ KA 

where , 

"Fg m Source or Cell Temperature 
T c ■ Mean Fluid Temperature 
Ax. « Distance, Fuel Cell — Interface 


°C 

°C 


m 
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4x 2 " Distance, Interface - Cooling Channel Wall o 

Observe that the value 1/hA is used to represent the 
convection resistance. The overall heat transfer by 
♦combined conduction and convection is expressed in terms of 
an overall heat-transfer coefficient U , defined by the 
relation , 


UA AT 


overall 


( 20 ) 


where A is the total heat transfer area. In accordance 

with Eq. (19) the overall heat transfer coefficient would 
be : 


U 


ITT 


i 


Ax ] 


( 21 ) 


K 


where , 

r c * thermal Contact Resistance °C.m/w 

Ax - Sum of Ax j and ax 2 m 

^ * Ciraphite Thermal Conductivity w/m.°C 

Thermal coatact resistance between the cooling plate and 
the fuel-ce-n plate was measured for various interface 
pressures and a correlation was developed relating the 
contact resistance with interface pressure. More details 
pretainii* to this correlation will be presented in an 
upcoming 5S^r.tion. 

5 • 5 Expeni tiu ntal Correlations 

MabhceMfcical correlations for the local heat transfer 

of POOfcHjMLRff 
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coefficient in terms of Nusselt number and Reynolds number 
for all coolants oil, water and air are presented. Also, 
the thermal contact resistance and the temperature drop 
^across the interface as a function of interface pressure 
are also shown . 

This arrangements permitted the presentation of the 
experimental results in their most general form in order to 
establish the basic relationships between the essential 
parameters, and to specify the significant parameters 
themselves . 

In an attempt to correlate the results (Appendix B), 
a computer program was developed to handle various types of 
curve fitting procedures including linear, polynomial, 
exponential, power and reciprocal functions. The least 
squares techniques and Gauss-Siedel approach, as well as 
other algorithms, were utilized in the development of this 
computer code. It also included graphics subroutine which 
enabled comparison of the various types of functions 
simultaneously on the screen of the graphic terminal. 

Thus, the experimental data shown in Tables 1.1 through 
1.5 were correlated by using the following functional 
variation , 


h - C Re m ?T n (22) 

where C, m and n are constants to be determined from 
experimental data. Figure 21 shows the dependence of the 
heat transfer coefficient on Reynolds number. By using all 


mmwL fm is 


Overall heal transfer Coeflicient, U c , w/m 






Figure 22- Variation, of the 
with Re and P fo 


oil as coolant 






Overall heal transfer Coefficient. U c , w/m 




the experimental data (Append ix B) as input to the computer 
program, a final value for m was determined as well as 
value for the constants n and C . 

Final correlations for all three coolants used -are 

<*■ 

represented as follows: 

(a) Serpentine Configuration (Pin 

The correlation equation which represents the heat 
transfer coefficient when oil is used as a coolant in terms 
of Nusselt number and Reynolds number is represented by, 

h - 16.444 Re °- 2144 Pr 0. 00769 (2J) 

and, t 


NU - 0.7501 Re 0 ' 2142 Pr °- 003 9 4 

( b > Serpen tine Configuration (Water) 

Also, similar equation was developed for the heat 
transfer coefficient when water is used as a coolant in 
terms of the same parameters as (a) , this equation is 
given by. 


h - 31.088 R e °- 0577 9 Pr 0. 01433 

• v - . ' ■ ' 


(25) 


aad , 


NU - 0.3259 R e ° • 05632 p r 0.0385 


(26) 


( c ) Straight Channel Configuration (Air) 

For the straight channel configuration and the cooling 
air, the heat transfer coefficient in terms of Nusselt 
number and Reynolds number can be written as follow^ 

h - 0.004338 Re 1 * 0730 p r 0 * 04085 ( 2 7) 


and , 


NU 


0.0005138 Re 


1.0966 


_ 0.237 
Pr 


( 28 ) 


(d ) Correlation of the Thermal Contact Resistance 

Similarly, correlations of the experimental data for 
thermal contact resistance, effective temperature drop 
across tke interface as a function of interface pressure 
were developed. Measured values of thermal contact 
resistanca versus interface pressure were plotted for the 
serpentine and the straight channel configurations as shown 
in Figure* 24 and 25 . These plots were then utilized to 

obtain tie following correlations: 

1- Straight Channel Configuration 

The following equation was developed for the contact 
resistance in terms of the interface pressure from the 
*_experimen.tel results. 


r 

e 


0.001581 


(-0.001413 P) 
e 


(29) 


is 

of pcaygmutt 


689 


1373 


2068 


INTERFACE PRESSURE . I pa 


Figure 24 - Variation of Contact Resistance with Interface Pressure for 
straight configuration 
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interface pressure, Kpa 


- Variation of Contact Resistance with Interface Pressure 
for serpentine configuration 



imteiface pressure f 


Figure 26 - Effect!?* Temperature Drop across the Interface versus 


Contact Pressure for serpentine and straight configurations 
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2- Serpentine Configuration 

Thermal contact resistance in terms of the interface 

»■ 

pressure was formulated for this configuration and at can 
^ written as follows: * 




0.00357 


(-0.002356 P) 
e 


(30) 


During this course of analysis and testings all the 
measurements for the temperatures used for the evaluation 
of the thermal contact resistance were taken without any 
coolant flowing through the cooling system. 

Also, experimental data for effective temperature drop 
for both serpentine and straight configurations were 
combined and plotted against interface pressure as shown in 
Figure 26 . The plot was then used to obtain the following 
correlation : 


AT - 7,986 - 0,9597 x 10 p kH) 

Results obtained from using the interval estimation 
technique (38) showed that these correlation equations have 
a confidence level of 95 percent. 

5.6 Thermal Contact Resistance Meas urements 

Figure 27 shows the locations of thermocouples that 
have ieen used for collecting temperature measurements for 
the thermal contact resistance on the cell plate and the 
cooling- plate. 


CWl PAfltJS 
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1 - Cooling Plate 

2 - Interface ni ^ 

3 - Fuel-Cell Plate 

4 - Thermocouples 

Figure 27 -Locations of Thermocouples at the Interface. 

Thermocouples located at positions shown in Figure 28 
were used to measure the temperature in the fuel-cell plate 
* during the test period and Figure 29 indicates the locations 
thermocouples used in the cooling plate. 

The determination of the temperature drop across the 
cell and the cooling plates interface was carried out by \*w. 

means of recording temperature measurements of thermocouples. 








contact between the two surfaces could be maintained. In 
ord er to ensure good contact between the cooling plate and 
the fuel-cell plate , graphite powder was injected to the 
contact area so as to minimize the existence of air 
initially contained in the unit. Application of pressure 
was implemented by the use of a hydraulic jack as shown in 
Figure 9 . The heat was determined from current and ' 

voltage measurements and the temperature drop across the 
interface was found by linear extrapolation as shown in 
Figure 30 . Thermal contact resistance was then evaluated 

by using Eq . (4) 

5 . 7 Analysis of Experim ental Results. 

In this section, an evaluation of the performance of 
the two types of cooling plates is presented. The 
variation of the overall heat transfer coefficient with the 
clamping pressure and Reynolds number for all three coolants 
is shown in Figures 34 through 37 . Also, Tables 1.1 to 

3.5 (Appendix B) were used to plot Figures 31 through 33 
for Nusselt number versus Reynolds number for oil, wator 
and air and for straight and serpentine configurations. 

1 - The Nu-Re Relationship 

The Nu-Re relationship for all coolants and the two 
types of cooling plates is shown in Figure 31 . The values 
of Nu varied from 4 to 8 for air, 1 to 3 for oil, and 0.4 
to 0.6 for water (Appendix B). It can be noted that air 
has produced the highest Nusselt number in comparison to 
6 il' -and water. This has been so because the straight 
Configuration has a large heat transfer area and the 
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2- The U-Re Relationship 

Tables 1.1 through 3.5 (Appendix B) were used to plot 
Figures 35 and 37 where the variation of the overall heat 
transfer coefficient with the clamping pressure and the 


Reynolds 

number 

(Re) is shown. 

It can 

be seen f 

rom 


Figure 35 

that 

the overall heat 

transfer 

coef f icie 

nt ha 

s 

increased 

as Re 

ynolds number and 

the clamping pres 

sure 


increased 

. Thi 

s can be clearly 

seen in 

Figures 35 

and 

37 

where P ■ 

0 Kpa 

and P « 3448 Kpa . 

The va 

lues of th 

e ove 

rail 

heat tran 
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oefficient for oil, water 
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about 20 

Z when 

the clamping pressure va 
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3448 Kpa . 

Water 

produced higher 
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heat trans 
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coefficient tha 

n oil and air. This was 

expected d 

ue to 

the 

nature of 

water 

for being a good 

coolant 

. Air at 
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ing 

pressare 

of 3448 Kpa gave higher 

value fo 

r the over 

all h 

eat 

transfer 

coef f i 

cient than water 

at Reyno 
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Ids number 

s ove 

r 

6000 . 
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In this se 

ction an overall 

evaluati 

on for all 

the 


experimental ob 

servations in rel 

ation to 

the cooli 

ng sy 

stem 

and tie two coo 

ling plates used 

in this 

experiment 

is 



implemented . 

In an attempt to evaluate each cooling plate 
configuration, it may be concluded that many factors should 
be considered in order to reach an optimum design for best 
coolitrg system. Such factors would be construction 
sinplicity, reliability, type of coolant and cost. The 
stirax>$bt channel configuration which includes the air 


of. foorwuiv 


cooling is by far the simplest in construction. Besides, 
in real fuel cell application, it keeps the cathodic 
reactant air separate from the cooling system gas (air). 
This improves electrical performance due to higher oxygen 
concentration of the cathode as well as reducing the need 
for acid resistant heat exchangers in the cooling 
stream (48) . 

The disadvantage of using the serpentine cooling 
system lies in its construction complexity and cost. It 
involves passing the liquid coolant through cooling plates 
inserted at regular intervals in the stack. The coolant 
may remain in the liquid phase at all times using only its 
sensible heat to cool the stack, or the coolant may 
partially vaporize using latent heat for part of the 
cooling load. The coolant must be either separately 
manifolded or supply lines must be connected to each 
individual cooling plate. Standard heat transfer materials, 
such as copper, are corroded by phosphoric acid and must be 
protected if it is to be used. Finally, the liquid itself 
and coolant lines, if electrically conductive, provide a 
possible shorting path in the cell (Shunt currents which 
result in parasitic losses) which must be minimized. 

However, in air cooling most air leaks could be 
ignored, but a leak in the stack liquid cooling system 
could cause a shutdown while the leak was repaired or the 
defective cooling plate replaced. 



Figure 3i -Nusselt Number 
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Figure 32 -Nusselt Number versus Reynolds Number for Straight and 
Serpentine Configurations with Air and Oil as Coolants. 
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CHAPTER VI 

EFFECTS OF HEAT TRANSFER COEFFICIENT 
ON ELECTRODE TEMPERATURE DISTRIBUTION 
AND OTHER STACK PARAMETERS 

The quantity of electrical energy produced in 
phosphoric acid fuel cells is accompanied by approximately 
equal amounts of generated heat energy. Removal of this 
excess heal energy is accomplished by the flow of reactant 
gases and the system of cooling plates. 

The feel cell heat generation is not evenly distributed 
across the cell plate because of depletion of reactant 
gases. This results in non-uniform temperature 
distribution. The experimental heating elements do not 
provide exact simulation of fuel cell heat generation, but 
they do load the forward edge of the plate at the coolant 

entrance. 

The work reported in this chapter is directed towards 
comparing the steady state temperature profiles for the 
experimental simulation of the fuel cell with the results 
obtained by Alkasab and Lu (6) in their computer simulatio 
model. K&o, effects of other parameters such as stack 


Serpentine Configuration 


A 
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pressures and coolant flow rates with air, water and oil. 

A mathematical correlation was developed to relate the 
average fuel-cell plate temperature with the air cooling 
flow rate and the stack clamping pressure. This correlation 
is represented by the following equation. 


a v 


max 


0.8916 - 0.0926 


P . m 

^max ' “max 


(32) 


where , 


T 

av 

Average cell 

T * 

A max 

190 °C 

^max" 

500 psi 

^min” 

0 psi 

“max* 

46.62 Kg/hr 

m min* 

0 Kg/hr 


Figure 38 shows the variation of the fuel cell plate average 
temperature as function of the above parameters. 

6.2 Effects of Cooling System Parameters on the 

Temperature Distribution 

The .temperature contour lines are very sensitive to 
many parameters such as the heating rate, coolant flow 
rate, stack clamping pressure and geometry of the cooling 
plate. 

Figures 39 *nd 40 show the effect of flow rate on 
temperature distribution when a serpentine cooling plate 
is used. Temperature measurements for the surface of the 


Straiphi Channels. Ajr 
i p - 0 Kpa 
» p - 689 kpa 
» p • 1379 kpa 
a p * 2068 kpa 
o p *2758 kpa 
o p *3448 kpa 


L experimental correlation 


P * m/puAl * m 

Figure 38 -Variation of Av.r.ge Fuel Cell Temperature with 
Coolant Flow Hate end Interface Pressure for 
Straight Configuration with Air as Coolant. 
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clamping pressure and coolant flow rate on the temperature 
distributions as well as the variations of the operating 
fuel cell temperature and system efficiency with the stack 
clamping pressure are discussed. 

6 . i Temperature Distributions 

In actual operation, the heat dissipation across the 
fuel cell plate is of asymmetrical type. Since only one 
curved boundary (the fuel cell plate) transfers energy to 
the flow. The experimental module was built in a manner 
where a real simulation of the actual operation was 
implemented . 

The temperature variations are an indication of the 
local heat generation rate where these are changed by the 
coolant flow rate. The temperature distribution for the 
fuel cell module was obtained by measuring the temperatures 
at seventy-two locations on the surface of the fuel-cell 
plate. Isotherms defining the temperature distribution 
across the surface of the fuel cell plate were drawn for 
various test cases including all three of the coolants used 
in this investigation. The average fuel cell temperature 
was calculated by obtaining the mathematical average of the 
temperature for all of the isotherms. 

These temperature profiles were a function of stack 
clamping pressure as well as the coolant flow rate. The 
clamping pressure showed to have a significant effects on 
the temperature distributions. 

Figures 69 through 74 (Appendix A) show the variations 
in the shapes of the isotherms for various stack clamping 
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Obviously, the increase of temperature is due to the 
influence of the clamping pressure on the system which 
resulted in producing a better contact between the cell 
plate and the cooling plate. Consequently, this resulted 
in a better temperature uniformity and a smaller temperature 
differential in the water cooling flow direction. Also, 
the average cell temperature increased from 168 °C to 174 °C 
which is within 2 percent of the average cell temperature 
obtained by Westinghouse measured data . 

Similar variations of the temperature distributions 
under the influence of stack clamping pressure and coolant 
mass flow rate have been also observed for oil and air. It 
can be noticed from Figures 41 and 42 that, when the stack 
clamping pressure increased from zero psi to 500 psi under 
the same coolant flow rate, the average cell plate 
temperature increased by 2 °C. Also, as the coolant flow 
rate increased, the temperature decreased by 15 °C at a 
flow rite of 88.20 Kg/hr. 

Another important result is that the contour lines of 
temperature corresponded to break closer to each other 
under the effect of higher stack clamping pressure. 

It may be concluded that most of the heat generated is 
removed through the cooling plate. But from a cost view 
point r the greater the flow rate of coolant the greater the 
auxiliary power needed to recycle the coolant and the more 
coolarc that is used, the lower the mean temperature will 
be, bat the greater the temperature difference between 
platesi will be in stack cooling situations. Lower peak 
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temperature meahs a more uniform temperature distribution 
and lover average temperature means that less auxiliary 
power will be required to pump the coolant. 

6.3 Thermal Analysis 

The energy balance equations for the fuel cell plate, 
cooling plate, process air and coolant were developed by 
Alkasab and Lv (6). In their analysis, they considered a 
phosphoric-acf < fuel cell stack similar to that shown in 
Figure 45 . It 1 such a stack, it was assumed that the 
number of fuel rie 11 modules contained between each pair of 
cooling plates to be five, and the repeating stack 
components to clhsist of a half cooling plate and two and a 
half fuel cell ^xfules. 



Figure 45 ^iaaet-rj of a Strip of Element for the Thermal 
t^B^dLs Model 
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The mathematical model developed by reference (6) 
includes four energy balance equations for the fuel cell 
plate, cooling plate, process air, and coolant are as 
follow. 

Fuel cell on process air side in air flow direction 

_ Cpmp 6Tp + (v* - V)I * 0 
x p p ay 


t Ky il + JL 

ty 1 ax 


Kx _6T 
x+t to 


Cooling plate in coolant direction 


t'Ky £ + 2 Kx £ 


ay 


6X 


x + t.’ 
2 


- C c m c 6l c » Q 
P c W 


Process air side 


dT 


* hp S (T 

dy mp C p 


V 


Coolant side 


- he Sc (T - T c ) 
dy me C c 

Boundary conditions 
x ■ 0 . »T/ ax - 0 

y - 0 aT/ ay - 0 

x - Lx dT/ to - 0 


(34) 


(35) 


(36) 


(33) 


symmetric condition 
adiabatic assumption 
symmetric condition 
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y . Ly 6T/ 6y « 0 adiabatic assumption 

y - 0 T p - T p , inlet 

y « 0 T c - T c , inlet 

where m « mass flow rate, Kg/hr-channel 
C - heat capacity, J/Kg-K 

Ky« effective thermal conductivity of cell in flow 
direction, J/hr-m-K 

K x « effective thermal conductivity of cell in 
stacking direction, J/hr-m-K 
t ■ thickness of cell including fuel and air channel, 
m 

x 1 m effective conduction distance from plate to upper 
cell plate, m 

x2 ■ effective conduction distance from plate to lower 
cell plate, m 
p ■ pitch of channel, m 

x 1 ’ * effective conduction distance from cooling plate 
to upper cell plate, m 

Lx, Ly - height and length of one slice, respectively, 

m 

V* - H/ZF , V 

t ' ■ thickness of cooling plate, m 
h ■ heat transfer coefficient, J/hr-m*-K 
S m perimeter of the channel, m 
Subscription 

p - process air 


c » coolant 


These simultaneous ordinary differential equations and the 
corresponding boundary conditions were solved by the finite- 
difference method. 

6>4 gp^tion of Experim ental Results-Computer Mod el 

The mathematical model of fuel cells that was 
described in the previous subsection, was developed into a 
computer program by reference (6) to facilitate performing 

the necessary calculations. 

This computer program was used in predicting the 
distribution of the temperature profiles in the fuel-cell 
stack. Experimental correlations developed, by the present 
investigation, and the overall heat transfer coefficient 
were incorporated into the computer program of reference (6) 
The modified program was then utilized to predict the 
temperature distribution in the fuel cell stack, and the 
results were compared with those obtained by reference (6). 

In considering the overall analysis for the cooling 
system, reference (6) considered the cooling factor to be 
a function of the heat transfer characteristics, plate size, 
and stack construction. The latter is primarily to specify 
the number of power plates between a pair of cooling plates. 
The he^t transfer characteristics is a function of the type 
of coolant (gas or liquid), cooling plate design, and the 
thermal conductivities of the plate material. In comparing 
the current experimental results with those available in 
the literature and reference (6), only the straight channel 
and the serpentine configurations were considered. 

- The developed computer model was modified to meet the 



cooling configurations used in this experiment. Also, 
factors such as thermal conductivities of plate material, 
plate size, and current densities were kept constant 
throughout the analysis. However, coolant mass flow rate, 
inlet and outlet temperature of the coolant, thermophysical 
properties and dimensions of the cooling channel were 
supplied to the program as input data. 

The program was run on the IBM 370 in NASA Lewis 
Research Center. The temperature of each grid on the cell 
plate, from outer-most plate to central plate was obtained. 
Also the average operating temperature, the operating 
pressure, and the DC output voltage for the specified stack 
were calculated. A difference between the results generated 
using the reference (6) data and the current experimental 
data is apparent. The difference is in the range of 5 to 
11 percent. 

6 . 5 Comparison of Experimental Results 

In this section, examination of the results obtained 
by the computer analysis describing major effects on the 
steady state temperature distribution and other aspects in 
relation to the effects of coolant and stack clamping 
pressure on the variations of temperature are presented. 

The computer outputs representing the temperature 
distribution on the fuel cell plate under the influence 
of the above parameters may le discussed as follows, these 
are representative outputs fr'un: 
a) Water Cooling 

. Figure 46 and 47 show the numeric distribution of tb« 
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Figure 46 -Computer Output for Temperature Distributions Generated by 
Experimental Data (Water Flow Rate ■ 47.50 Kg/hr , Current 
Density 0.325 am/cm* , Stack Clamping Pressure « 0 Kpa) 
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1 THE AVERAGE OPERATING TEMPERATURE IS 8.44543E 03 K 
THE OPERATING PRESSURE IS J.O ATM 
THE PULL DC POWER OUTLET IS 8.52974E 11 KU-DC 


Figure 47 -Computer Output for Temperature Distributions Generated by 
Experimental Data (Water Flow Rate « 47.50 Kg/hr , Current 
Density 0.325 am/cm 1 , Stack Clamping Pressure ■ 3448 Kpa 


temperature in each finite difference on the cell plate at 
the final steady state under the Influence of the stack 
clamping pressure. 

.Figures 46 and 47 were generated by the experimentally 

A -.. 'M ,* r/;f 
"V " - • 
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a .,. lop .d heat transfer coefficient, with «t.r cliM »< 
t t ero end 3448 Kpa .tact clamping pressures. The eff.cta 
*, th , .tack clamping pressure la clearly noticed since the 
average fa.l cell pl»» temperature decreased from 
„ l76 -C when the stack clamping pressure increased from 


zero to 3448 Kpa. 
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Figure IS -Computer 

SSToSs W«F Tstack Clamping Presaure - -3448 Kpa, 
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Comparison between reference (6) data and the current 
experimental data is shown in Figures 48 and 49 , these 
temperatures were produced using water as coolant and with 
equal flow rates. The outputs of these temperatures showed 
an average cell temperature of 186.10 °C for the 
experimental data, while the reference (6) data for the 


same flow rate yielded an average fuel cell temperature of 
204 °C. A difference of 11 percent is apparent between 
experimental and predicted results by (6). 


b ) Oil Cooling 

Similarly for oil cooling, Figures 50 and 51 represent 
the numeric outputs of temperatures at the final steady 
state . 
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THE OPERATING PRESSURE IS l.o ATM 

THE FULL 0C POWER OUTLET IS 0.59286E 01 MJ-DC 


Figure 50 -Computer Output for Temperature Distributions Generated by 
Experimental Data (Oil Flow Rate * 16.66 Kg/hr , Current 
Density 0.325 am/cm* , Stack Clamping Pressure - 0 Kpa) 


The temperature of the fuel cell plate decreased from 
200 °C to 193 °C when the stack clamping pressure Increased 
front zero to 3448 Kpa. Also, the average fuel cell plate 
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temperature varied between 190 °C and 200 °C at various 
coolant flow rate for the experimental data and between 
175 °C and 204 °C for reference (6) data. The differ nee 
is about 10 percent. 
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Figure 51 -Computer Output for Temperature Distributions Generated by 
Experimental Data (Oil Flow Rate - 16.66 Kg/hr , Current 
Density 0.325 am/cm J , Stack Clamping Pressure -3448 Kpa) 
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1 THE AVERAGE OPERATING TEMPERATURE ZS 0.46924E 03 K 
THE OPERATING PRESSURE IS 1 ATM 
THE FULL DC POWER OUTLET IS 0.58043E 01 KU-DC 


Figure 52 -Computer Output for Temperature Distributions Generated by 
Experimental Data (Oil Flow Rate - 24.87 Ig/hr , Current 
Density 0.325 am/ cm* , Stack Clamping Pressure ■ OKpa) 
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1 THE AV«a CE OPERATING TEMPERATURE IS 0.44579E 03 K 
THE OPERATING PRESSURE IS 1 ATM 
THE FULL K POUER OUTLET IS 0.53028E 01 KU-DC 


Figure 53 -Co«|uter Output for Temperature Distributions Generated by 
Refeience (6) Data (Oil Flow Rate - 24.87 Kg/hr , Current 
Density 0.325 am/cm* , Stack Clamping Pressure - 0 Kpa) 

c ) Air Cooling 

For air cooling steady state temperature distributions 
are r epr e seated by Figures 54 and 55 under zero and 500 psi 
stack clamping, pressures. 

The avenge fuel cell plate temperature dropped from 
192 °C to 17-6 *'C as the stack clamping pressure increased 
from zero to 3448 -Kpa. Also, the temperature distributions 
showed some difference under the influence of the heat 
transfer cteff-icients between experiment and predicted by 
(6). Neveckaless , the agreement was within 2 percent 
which is an-' lid i cat ion of validity for the computer 
simulation The average fuel cell temperature was about 
190 °C fortatth experimental and reference (6) data. 
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Figure 54 —Computer Output for Temperature Distributions Generated by 
Experimental Data (Air Flow Rate - 29.84 Kg/hr , Current 
Density 0.325 am/cm’ , Stack Clamping Pressure - 0 Kpa) 
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Figure 55 -Computer Output for Temperature Distributions Generated by 
Experimental Data (Air Flow Rate -29.84 Kg/hr , Current 
Density 0.325 am/cm* , Stack Clamping Pressure - 3448 Kpa) 
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Figure 56 -Computer Output for Temperature Distributions 
Experimental Data (Air Flow Rate ■ 29.84 Kg/hr 
Density 0.325 am/cm* , Stack Clamping Pressure 
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Figure 57 -Computer Output for Temperature Distributions Generated by 
Reference (6) Data (Air Flow Rate ■ 29.84 Eg/hr , Current 
Density 0.325 am/cm* . Stack Clamping Pressure - O Kpa) 


The efficiency^, of the fuel-cell stack can be 
defined as, 

_ E s - E p (37) 

* E + Q 
*■ s y s 

where , 

E g ■ Eltctric energy generated by the fuel-cell stack 
Ep m Ele> trie energy consumed by the pump to circulate 
the coolant. 

Q s ■ Heat generated by the fuel-cell stack 


where E p is directly proportional to the coolant mass flow 
rate. On the ot'ier hand, it can be shown that the mass 
flow rate is inversely proportional to the value of the 
overall heat transfer coefficient U c . From equations (11) 
and (20) the following relation can be obtained, 

la. ~ T * . _L i ♦ _L._L 08) 

0 A P T c 2C P " 

where T g , , q, A p and C p are assumed to be constant. 

Eq. (38) can also be written as, 
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where , 



Values of m were calculated for a range of values of U c 
for air, oil and water, and the results were plotted against 
the experimental values of U c and their corre^oonding m's . 
These plots are shown in Figures 58 , 59 and 60 respectively. 

It can be noted from these figures that the maximum 
deviation between the calculated and the measured values of 
U c and m is about 14 percent for water, 16 percent for oil 
and 9 percent for air cooling. 

Also, these figures clearly indicates that the higher 
the value of U c , the lower is the value of the flow rate m. 
Since the work consumed by the pump to circulate the 
coolant is directly related to the flow rate, Eq. (37) 
shows that the higher is the value of U c , the greater will 
be the efficiency of the fuel cell stack. Also, since U c 
is directly proportional to the clamping pressure, one can 
conclude that the higher is the clamping pressure, the 
higher is the efficiency of the fuel-cell stack. 
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Figure 58 -Coolant Flow Rate Requirements. To maintain The Electrodes 
Temperature at 190 °c, As a Function Of The Overall Heat 
Transfer Coefficient . (Serpentine Configuration, Water Cooling) 
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Figure 59 -Coolant Flow Rate Requirements, To Maintain The Electrodes 
Temperature at 190 °c, As Function Of The Overall Heat 
Transfer Coefficient. (Serpentine Configuration, Oil Cooling) 
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Figure 60-Coolant Flow Rat, Retirements, To maintain The Electrode, 
Temperature at 190 •c.Aa a function Of The Overall Heat 
Transfer Coeficient. (Straight Configuration, Air Cooling) 




CHAPTER VII 

RESULTS AND DISCUSSIONS 

Some of the important features pertaining to the 
results obtained during this course of investigation are 
discussed below. 

While the investigation was primarily concerned with 
heat transfer, preliminary measurements of the thermal 
contact resistance were also made to evaluate its influence 
on the overall heat transfer coefficient. The values of 
the overall heat transfer coefficients for the three 
different coolants and the two different cooling plate 
configurations are plotted in Figures 6l , 6l and 63 . In 
all cases, the overall heat transfer coefficient increased 
with increasing interface pressure as shown in the figures. 

The interface pressure was varied from 0 Kpa to 3448 Kpa in 
incremears of 689 Kpa . For oil cooling in Figure 61 the 
flow wasin the laminar range. Between Re - 22 and Re - 43 
the rateof increase of U c is doubled due to the doubling 
of theflVw rate. Figure 62 shows the variation of the 
overalU hemt transfer coefficient for water cooling. The 
transient. of flow from laminar to turbulent occured 
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between Re - 1740 and Re - 3321 . In Figure 63 for air 
cooling, the transition from laminar to turbulent is 
between 3070 and 4341. 

Figure 64 show the relationship of the ratio of the 
overall heat transfer coefficients without and with the 
thermal contact resistance term (U/U c ) decreasing as the 
pressure increases. (U/U c ) varied from 1.15 at P ■ 0 Kpa 
to 1.05 at P « 3448 Kpa , for oil and water, while for air it 
varied from 1.05 to 0.7 at P * 0 Kpa and P ■ 3448 Kpa 
respectively . 

Two series of test were conducted for measurements of 
the thermal contact resistance and the data obtained was 
plotted for both configurations, serpentine and straight. 

As shown in Figure 25 the thermal contact resistance varies 
from 0.0016 at 0 Kpa to 0.000812 m*.°C/w at 3448 Kpa . 

Figure 65 and 66 show the effect of increasing interface 
pressure on the thermal contact conductance. As expected, 
the conductance increases with increasing contact pressure. 
The values of contact conductance measured during the 
unloading phase are significantly higher than the values 
obtained during loading phase. This increase of conductance 
can be explained by the better conformity of two interfaces 
after the maximum load has been reached, i.e., a plastic 
deformation of the surface asperities has resulted in a 
greater true area of contact. 

The effective temperature drop across the interface of 
the cell/cooling plates for both cooling plate 
configurations was plotted versus pressure, in Figure 26 , 


\T 


INTERFACE PRESSURE , P , K P a 

Figure 64 - U/U c versus Interface Pressure for Serpentine and 
Straight Configurations with Air, Water and Oil as 


Coolants. 




straight configuration 
• last no. 1 
▼ last no. 2 



Figure 66- Effect* of Contact Pressure on the Thermal Conductance 


for straight configuration 


tb . numerical correlation for this data is represented b, 

' (31) . Similarly. Eq. (29) (30) represent the 

contact resistance as function of pressure for the two 

cooling plate configurations. 

The experimental heat transfer condition is one of 
•symmetrical heating since only the curved boundary (i.e. 
cell plate), transfers energy to the flow. If the average 
experimental heat transfer coefficient for this condition 
is used to compare the corresponding experimental 
temperature data, the nature of the thermal boundary 
condition is irrelevant. Nusselt number values for all 
three situations were plotted versus Reynolds number on a 
log-log scale. For each case, data for three tests were 

plotted as shown in Figures 67 , 68 and 69 . 

The data for the heat transfer measurements such as 
heat transfer coefficient and Nusselt number shown in 
Tables 1.1 through 3.6 were correlated into a single 
relation, which is given by Eq . (22). Consequently, 

equations representing data for oil. water and air cooling 
were formed and results obtained by these equations were 
plotted against the experimental results as shown in 
Figures 67 , 68 and 69 . Table 2 represents measured 
measured and calculated data by equations (24), (26) and 
(28) for all three coolants, oil, water and air. 

The experimental error could be reduced by maintaining 
the largest possible temperature difference between the 
coolant and the cell plate. Since the heat transfer rates 
increased with increasing Reynolds numbers, the input power 


to the cell plate had to be limited to keep the cell 
temperatures below 200 °C to prevent damage to the 
thermocouple adhesive bonds. This condition was used for 
B ost of the tests. There was no significant change in the 
heat transfer rate when the power dissipated in the cell 
was reduced by one half. This showed that the inaccuracies 
in the thermocouple measurements did not significantly 
effect the results. 


TABLE 2 

Comparison Between Experiment and Correlation 
Air Case 


Re 

3070 

4341 

5317 

6139 

6864 

NU exp 

3.005 

4.840 

6.872 

6.632 

7.28 

NU C0 

3.148 

4.61 

5.744 

6.72.4 

7.60 

Water Case 

Re 

1165 

1740 

3321 

4587 

6167 

NU 

0.525 

0.520 

0.545 

0.55 

0.565 

exp 

NU 

0.488 

0.499 

0.517 

0.527 

0.545 

CO 

Oil Case 

Re 

15 

22 

43 

59 

79 

NU 

1.34 

1.49 

1.74 

1.86 

1.93 

exp 

NU co 

1.373 

1.491 

1.723 

1.842 

1.961 


where , 

NU - Nusselt Number (Experimental) 
exp 

NU , ■ Nusselt Number (correlation) 

CO 
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Thr « measurements for the h..t transfer co.fflci.nt. 
re .ad« for each coolant under the .ana conditions. No 
’ ‘ differences were .ad, for each coolant under the same 

“adltions. No sajor differences -ere noticed except that 
for each coolant a different heat transfer rat. 
vhich ... expected. The greatest heat transfer rate -as 

obtained .ith air cooling. The air exit te.perature 
averaged about 60 -C -hil. for oil and .star cooling it 
averaged about 120 and 60 -C respectively. The Neynolds 
number ranged from 15 to 80 for oil. For air and water 
ranged fro. 1000 to 7000 which is -ell in the turbulent 

region for air and both la.inar and turbulent for water. 

r _ _ rate was obtained with water 
The lowest heat transfer r 

cooling. Comparing the oil and water cooling cases, 
noticed that the values of the heat transfer coefficient 

were nearly agual. -.pie ° £ 

these results, consider the design of a nuclear reactor 

. . „ i S a known function of 

cooling system where the heat flux is a X 

V 0 Tvnieally the heat flux is low near 
length along the tube. \Typical y 

the entrance and exit reaching a .axi.um at the midpoint 
If the coolant isaliguid metal, considerable error in the 
predicted surface temperatures can result unless van 
heat-flux theory is employed. On the other hand, if the 
coolant is a gas. or pressurised water, the varying heat 
flux has little influence, and it is ad.gu.t. to use • 
Nusselt number based on constant heat flux theory to 
calculate the local temperature difference between the 

fluid and the wall aurfaces. 


In general terms, the accuracy of the experimental 
data and the results, is considered good, although some 
factors such as surface roughness, friction, geometry of 
the cooling channels and configuration of plates could have 
effected the results to some extent. The effect of the 
surface roughness on the heat transfer for flow in tube is 
essentially the same as that for the external boundary 
layer. The only differences arise from the fact that for 
flow in tube the friction coefficient is based on the mean 
flow velocity and the heat transfer coefficient is based on 
the mixed mean fluid temperature. 

It should also be noted that artificial roughness or 
turbulence promoters are frequently employed as a means to 
increase the heat transfer coefficient. A comprehensive 
review of methods for augmentation of convective heat 
transfer is given by Bergles (35), (36) • However, 
differences between the measured average heat transfer 
coefficient and the predicted by the use of existing 
correlations do exist. This difference is a result of the 
effects of asymmetry of heat transfer and the inadequacy of 
the hydraulic radius concept itself in the case of 
noncircular cooling channels. The average heat transfer 
coefficient used in the theoretical prediction was taken 
from existing correlations, and of course that correlation 
refered to a particular geometry and thermal boundary 


conditions. 

A tabulation of measured heat transfer coefficients, 
thermal contact resistance and other parameters are all 
presented in the appendix. 
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Figure 68 -Naselt Number versus Reynolds Number for Straight 
Configuration with Air as a Coolant. 
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ure 69 -Nusselt Number versus Reynolds Number for Serpentine 
Configuration with Oil as a Coolant. 
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CHAPTER VIII 

CONCLUSIONS AND RECOMMENDATIONS 
M CONCLUSIONS 

A. Based on the experimental results, it may be con, 
that each of the cooling system used in this experiment he. 
its advantages and disadvantages as far as the performance 
of the cooling system is concerned. Performance 
characteristics for each of the cooling plate can be listed 

as follows: 

1. The advantages of using the straight channel configuration 
lies in its special characteristics of having a short flo 
length, large effective heat transfer area, good temperature 
uniformity and smaller temperature differential between 
coolant and cell plate. Also, it has a high reliability and 


ow cost. 

However, some disadvantages are summerized by its 
imitations to the use of gases only which limits the size 
,f the plate since gases have low heat capacity. Besides the 
lif f iculties in controling air leakage. 

!. The serpentine configuration does not have any limitations 
towtds the size of the plate (i.e. plate can be large) since 
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liquids have high heat capacity. Another advantage is 
the capability of using two phase flow coolant (41). 

Nevertheless, the serpentine configuration is costly and 
complex in construction. Also, it provides less temperature 
uniformity, high inlet and outlet coolant temperature 
differential. It has a relatively small heat transfer area, 
long flow path which makes it not applicable for gases. 

B. the cooling system performance and the efficiency of 
the fuel-cell stack is a strong function of the clamping 
pressure. The effect of the clamping pressure for the two 
configurations investigated can be summerized as follows: 

1. For both configurations straight and serpentine with 
water, oil and air cooling, the higher the clamping pressure, 
the higher the value of the overall heat transfer coefficient, 
the greater will be the efficiency of the fuel-cell stack. 

2. Higher heat transfer coefficient has been obtained for 
air cooling, while it has been nearly equal for oil and 
water cooling. 

8 . 2 Recommendations 

The present study is the first in correlating heat 
transfer coefficient for the cooling system in the fuel-cells. 
Future studies are recommended in the following directions. 

1. Studies of the transient effects of the cooling system 
parameters on the overall heat transfer coefficient and the 
uniformity of electrode plate temperature distribution* 

2. Development of a mathematical relation for the coolant 
flow rate as a function of electric load, the stack clamping 
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TABLE 3 

•Comparison of the Analytical and Experimental Results (Water Case) 


Re 


T 

out 

(Measured ) 

T 

out 

(Calculated) 

P - 0 

psi 

(0 Kpa) 


1165 


48.40 

53.710 

1740 


47.50 

52.366 

3321 


46.00 

47.598 

4587 


45.20 

46.560 

6165 


44.40 

43.923 

P » 100 psi 

(689 Kpa ) 


1165 


51.10 

54.567 

1740 


49.53 

51.278 

3321 


46.95 

47.775 

4587 


44.75 

44.732 

6165 


43.10 

43.553 

P - 200 

psi 

(1379 Kpa) 


1165 


55.20 

55.894 

1740 


53.76 

52.248 

3321 


48.74 

47.735 

4587 


44.31 

* 

45.138 

6165 


46.29 

44.902 
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TABLE 3 - Continued 

♦Comparison of the Analytical and Experimental Results (Water Case) 


Re 

T . 
out 

(Measured) 

T . 
out 

(Calculated) 

P - 300 Dsi 

(2068 Kpa ) 


1165 

54.20 

54.737 

1740 

52.41 

51.079 

3321 

45.35 

45.723 

4587 

45.79 

45.072 

6165 

43.10 

42.832 

P - 400 psi 

(2758 Kpa) 


1165 

55.40 

55.849 

1740 

55.54 

52.360 

3321 

51.76 

48.469 

4587 

53.50 

47.757 

6165 

55.44 

47.548 

P ■ 500 psi 

(3448 Kpa) 


1165 

47.32 

49.760 

1740 

45.52 

47.463 

3321 

39.31 

42.106 

4587 

38.51 

40.717 

6165 

39.93 

40.813 
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TABLE 4 

•Comparison of Analytical and Experimental Results (Air Case) 


p , Q psi (0 Kpa) 


Re 

T - 
out 

(Measured) 

T , 
out 

(Calculated) 

3070 

71.40 

73.649 

4341 

66.38 

68.131 

5317 

65.62 

64.537 

6139 

63.96 

58.272 

6864 

61.00 

56.791 

P ■ 100 psi 

(689 Kpa) 


3070 

73.40 

75.413 

4341 

66.38 

66.042 

5317 

58.62 

55.924 

6139 

63.96 

56.073 

6864 

53.50 

50.1974 

P - 200 psi 

(1379 Kpa) 


3070 

74.10 

70.031 

4341 

67.38 

64.497 

5317 

65.22 

59.943 

6139 

63.00 

54.692 

6864 

61.34 

51.061 


TABLE 4 - Continued 

•Caparison of Analytical and Experimental Results (Air Case) 


Re 

T 

out 

(Measured) 

T - 
out 

(Calculated) 

P - 300 psi 

(2068 Kpa) 

) 

3070 

73.20 

75.385 

4341 

70.40 

74.293 

5317 

66.78 

63.787 

6139 

65.20 

59.193 

6864 

64.55 

54.511 

P = 400 osi 

(2758 Kpa) 


3070 

73.87 

67.146 

4341 

67.57 

69.981 

5317 

63.04 

61.292 

6139 

55.24 

48.929 

6864 

56.67 

48.379 


P - 500 psi 

3070 

4341 

5317 

6139 


(3448 Kpa) 

77.00 

68.10 

64.00 
61.81 
63.47 


70.769 

67.467 

60.669 

56.460 

59.924 


6864 
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TABLE 5 

•Comparison of the Analytical and Experimental Results (Oil Case) 


Re 

T „ 
out 

(Measured ) 

T „ 
out 

(Calculated ) 

P - 

0 psi (0 Kpa) 


15 

107.50 

107.333 

22 

106.60 

107.369 

43 

90.70 

88.423 

59 

75.70 

77.644 

79 

65.20 

70.298 

P = 

100 psi (689 Koa) 


15 

108.20 

114.830 

22 

107.90 

105.649 

43 

92.30 

88.871 

59 

82.65 

79.678 

79 

66.85 

69.777 

P - 

200 psi (1379 Kpa) 


15 

107.80 

113.362 

22 

111.65 

106.820 

43 

95.60 

89.957 

59 

80.15 

77.983 

79 

67.20 

68.738 



g. 

15 

98.57 

107.369 

i 




f 

22 

102.80 

101.085 

* 

$ 

43 

91.00 

84.988 

- 

59 

78.30 

75.010 

t§ 

79 

75.75 

71.819 

"^r 

P = 400 

£si (2758 Kpa) 


j 

15 

95.10 

103.888 

5 

22 

101.75 

99.599 


43 

90.45 

83.245 


59 

78.40 

74.489 


79 

65.65 

66.053 



P » 500 psi (3448 Kpa) 


15 

90.42 

100.48 

22 

95.32 

94.942 

43 

85.25 

80.692 

59 

74.14 

71.295 

79 

62.639 

61.120 


TABLE 6 

•Contact Resistance Measurement for Straight Channels Cooling Plate 


Test #1 


Pressure 1 

(psi) , (Kpa) ( C) 


r . 
c * 


(°C.ra 1 /w) 


h - 1/r 
c 

(w/°C.m I 


100 

689 

6.8 

0.001499 

667.11 

200 

1379 

5.6 

0.001234 

810.37 

250 

1724 

5.0 

0.001102 

907.44 

300 

2068 

4.7 

0.001036 

965.25 

400 

2758 

4.0 

0.000882 

1133.79 

450 

3103 

3.62 

0.000805 

1246.88 

500 

3448 

3.51 

0.000780 

1282.05 



Test #2 



100 

689 

6.7 

0.001478 

676.59 

200 

1379 

5.4 

0.001195 

836.82 

250 

1724 

3.65 

0.000810 

1234.57 

300 

2068 

5.1 - 

0.001125 

888.89 

350 

2413 

3.9 

0.000860 

1162.79 

400 

21 SS 

4.44 

0.000980 

1020.41 

450 

3 ICE.’ 

3.8 

0.000840 

1190.48 

500 

344fr 

4.0 

0.000880 

1136.36 


TABLE 7 


•Contact Resistance Measurements for Serpentine Cooling Piece 

Test #1 


Pressure 

(psi), (Kpa) 
\ 


200 

250 

300 

400 

500 


1379 

1724 

2068 

2758 

3448 


T 

(°C) 


6.1 

5.3 
4.75 
3.55 

3.3 


(°C.m7w) 


100 

689 

8.0 

0.00281 

200 

1379 

6.5 

0.00228 

300 

2068 

4.8 

0.00168 

400 

2758 

4.0 

0.00140 

500 

3448 

3.1 

0.00109 



Test #2 


100 

689 

7.8 

0.00275 


0.00215 

0.00185 

0.00167 

0.00125 

0.00115 


h c - l/r c 

(vrcV 


355.87 

438.60 

595.24 

714.28 

917.43 


363.64 

465.12 

540.54 

598.80 

800.00 

869 . 56 



34.82 


36.66 


38.03 


1.59 


1.67 


30.62 


32.04 


88.20 


Test #1, Oil with P= 0 psi ( 0 kpa) 


1.73 33.08 


22 . 

28 

24 . 

38 

27 

.60 

28 

.75 

29 

. 58 


TABLE 1.1 -Continued 


(Kg/Hr) 

(°C) 

(°c) 

16.66 

129.97 

86.50 

24.87 

181.10 

85.60 

47.50 

151.53 

69.70 

65.60 

132.43 

54.70 

88.20 

119.13 

44.20 
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31 .45 


35.22 


38.17 


39.49 


Nu 

U 

( w/m^ * ° C ) 

1.29 

25.46 

1.43 

27.98 

1.60 

30.93 

1.74 

33.18 

1 .80 

34.18 



23.76 


25.94 


28 . 44 


30.34 


31.17 
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TABLE 1.3 



Re 

h 

( w/m 2 ■ °C ) 

Nu 

U 

( w/m 2 • °C ) 

Uc 

(w/m 2 .°C) 

1 6 . 66 ~ 

15 

28.96 

1 . 3 2' 

25.99 

24 . 57 

24.87" 

22 

32.03 

1.46 

28.44 

26.75 

47.50 

43 

35.76 

1.62 

31.35 

29 . 29 

i 

65.60 

59 

38.61 

1.76 

33.51 

31 .18 

88.20 

79 

40.86 

1.86 

35.19 

- 

32 . 64 


Test #1, Oil with P- 200 psi( 1379 kpa) 


TABLE 1.3 -Continued 


m 

(Kg/Hr) 

Tw 

(°C) 

Tc 

(°C) 

16.66 

185. 20 

86.80 

24.87 

179.61 

90.65 

47.50 

154.27 

74.60 

65.60 

132.95 

59.15 

t 

88.20 

115.93 

46.20 































TABLE 1.4 


(Kg/Hr) 


16.66 


24.87 


47 . 50 


65.60 


88.20 


h 

(w/m 2 -°C) 

Nu 

29.39 

1.34 

32.67 

1.49 

38.17 

1.74 

40.82 

1.86 

42.34 

1.93 


(w/m 2 -°C) (w/m 2 -°C) 


26.34 


28.95 


33.18 


35.17 


36.29 


Test # 1, Oil with P= 300 psi (2068 kpa) 


25.17 


27.54 


31.35 


33.12' 


34.11 


TABLE 1.4 -Continued 


m 

(Kg/Hr) 

16.66 


24.87 
47.50 
65.60 
88.20 


(°C) 


(°C) 


174.50 77.57 

168.97 81.80 
144.65 70.00 
127.10 57.30 
122.05 54.75 
























6.66 


( w/ m 2 * ° C ) 


30.27 


33.34 


39.67 


43.21 


Nu 


1.38 


1.52 


1.81 


1.89 


1.97 


u 


( w/m 2 

•°c) 

26. 

34 

28. 

95 

33. 

,18 

35. 

.17 

36 

.29 


Test #1. Oil with P* 400 psi (2758 kpa) 



33.12 


3 


TABLE 1.5 -Continued 


(Kg/Hr) 

(°C) 

(°C) 

16.66 

168.28 

74.10 

24.87 

166.20 

80.75 

47.50 

141.27 

69.45 

65.60 

126.07 

• 

57.40 

fift . 20 

110.58 

44.65 
























/ 
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(Kg/Hr) 


6.66 


24.87 


.50 




65.60 

59 

88.20 

79 


30.71 


34.22 


39.53 


43.00 


44.75 


Nu 

U 

(w/m 2 -°C) 

1.40 

27.39 

1.56 

30.16 

1.80 

34.21 

1.96 

36.77 

2.04 

38.05 


Test #1 , Oil with P= 500 psi (3448 kpa) 


TABLE 1.6 -Continued 



26.60 


29.19 


32.98 


35.36 


36.53 


m 

(Kg/Hr) 

T w 

(°C) 

Tc 

(°C) 

16.66 

162.20 

69.42 

24.87 

157.58 

74.32 

47.50 

136.33 

64.25 

65.60 

119.41 

53.14 

88.20 

103.79 

40.12 
























A 


TABLE 2.1 


( Kg/ Hr )- 


6.66 


24,87 


47.50 


65.60 


88.20 


1165 


1740 


3321 


4587 


6167 


44.74 | 0.462 


45.13 0.466 


46.87 0.484 


47.17 0.487 


51.33 0.530 


Test #1, Water with P* 0 psi (0 kpa) 
TABLE 2.1 -Continued 


38.04 


38.32 


39.57 


39.78 


42.70 


33.49 


3 


34.83 


m 

(Kg/Hr) 

Tw 

(°C) 

Tc 

(°C) 

16.66 

91 .08 

27.40 

24.87 

92.46 

26.50 

47.50 

85.77 

25.00 

65.60 

84.61 

24.20 


88.20 


78.91 


23.40 
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TABLE 2.2 


m 

(Kg/Hr) 



1165 


1740 


3321 


4587 


88.20 6167 


( w/m^ *°C) 


45. 

52 

46. 

49 

47 , 

.36 

50 

.65 

51 

.04 


0.470 


0.480 


0.439 


0.523 


0.527 


I ( w/m^ , 0 C ) ( w/m^ • °C ) 


38.60 


39.29 


39.92 


42.50 


Test #1, Water with P* 100 psi (689 kpa) 


TABLE 2.2 -Continued 



35.88 


37.72 



m 

(Kg/Hr) 


(°C) 


16.66 92.69 

24.87 89.82 
47.50 86.11 
65.60 80.00 
88.20 77.93 


Tc 

(°C) 

30.10 


28.53 

25.95 

23.75 

22.10 
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m 

(Kg/Hr) 

Re 

h 

(w/m 2 -°C) 

Nu 

U 

( w / m 2 • 0 C ) 

Uc 

( w/m 2 • 0 C ) 

16.66 

1165 

46.49 

0.480 

39.29 

36.13 

24.87 

1740 

48. 10 

0.496 

40.44 

37.09 

47 . 50 

3321 

48.91 

0.505 

41.01 


65.60 

4587 

49 . 39 

0.510 

41.35 

37.86 

88.20 

6167 

50.85 

0.525 

42.36 

38.71 


Test #1, Water with P* 200 psi (1379 kpa) 


TABLE 2.3 -Continued 


m 

(Kg/Hr) 

Tw 

(°C) 

Tc 

(°C) 

16.66 

95.49 

34.20 

24.87 

92.00 

32.76 

47.50 

86.00 

27.74 

65.60 

81.00 

23.31 

88.20 

81.33 

25.29 




































24.87 


1740 


49.20 


0.492 40.12 


0.508 41.22 


37.47 



65.60 


3321 


4587 


50.21 


50.85 


0.517 


0.525 


88.20 


6167 


52 . 76 


0.544 


Test #1, Water with P» 300 psi(2068 kpa) 


TABLE 2.4 -Continued 


m 

(Kg/Hr) 


(°C) 


<°C) 


16.66 93.00 


33.20 


24.87 89.33 


31.41 


47.50 81.10 


24.35 


65.60 80.83 


24.79 


88.20 76.10 
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16.66 


24.87 


47.50 


65.60 


88 . 20 



1165 


1740 


3321 


4587 


6167 


h 

( w/m 2 • °C ) 

Nu 

48.33 

0.499 

49.39 

0.510 

49.98 

0.516 

51 .81 

0.535 

53.20 

0.549 



40.60 


41.35 


41.76 


43.03 


43.99 


Test #1, Water with P* 400 psi (2758 kpa) 
TABLE 2.5 -Continued 


(w/m 2 -°C) 


38.43 


39.10 


39.45 


40.60 


41.45 



m 

(Kg/Hr) 


16.66 

24.87 

47.50 

65.60 


88.20 


(°C) 


95 . 36 
92.23 
87.77 

87.50 


88.00 


(°C) 


36.40 

34.54 

30.76 

32.50 


34.44 
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TABLE 2.6 


m 

(Kg/Hr) 

Re 

h 

( w/m 2 • °C ) 

Nu 

U 

(w/m 2 -°C) 

Uc 

(w/m 2 .°n 

16.66 

1165 

• 

50.85 

0.525 

42.37 

40.50 

24.87 

1740 

50.36 

0.520 

42.03 

40.19 

47.50 

1 

332 1 

52.78 

0.545 

43.70 

41.71 

65.60 

4587 

53.27 

0.550 

44.03 

42.02 

88.20 

6167 

54.72 

0.565 

45.02 

42.91 


Test #1, Water with P« 500 psi (3448 kpa) 


TABLE 2.6 -Continued 


m 

(Kg/Hr) 

Tw 

(°C) 

Tc 

(°C) 

16.66 

82.36 

26.32 

24.87 

81 . 10 

24.52 

47.50 

72.30 

18.31 

65.60 

70.00 

17.51 

88.20 

71.00 

18.93 
















TABLE 3.1 


m 

(Kg/Hr) 


0.6950 3070 


0.9829 4341 


1.2038 5317 


( w/m 2 *°C ) 


19.966 


32.518 


43.713 


1.3900 6139 . 43.99 


1.5540 I 6864 I 48.375 


3.005 


4.895 


6.5799 


6.622 


7.2817 


(w/m 2 -°C) (w/m 2 -°C) 


18.51 


28.83 


37.29 


37 . 50 


40.63 


Test #1, Air with P= 0 psi k P a ' 


17.98 


27.57 


35.2 


35.40 


38.18 


TABLE 3.1 -Continued 


m 

(g/Hr) 

Tc 

(°C) 

ATc 

(°C) 

ATm 

(°C) 

0.6950 

50.10 

55.00 

74.38 

0.9829 

48.41 

50.19 

58.94 

1.2038 

45.60 

45.20 

53.70 

1.3900 

41.22 

35.95 

44.13 


1.5540 


39.95 


33.70 


42.06 





































TABLE 3.2 


m 

(Kg/Hr) 


0.6950 3070 


0.9829 4341 


1.2038 5317 


1.3900 6139 


1.5540 | 686 


Test #1 , Air 




20.156 


34.500 


41.08 


.49.48 


49.175 


3.034 


5.193 


6.1836 


7.448 


7.402 


u 

( w/m 2 *°C ) 

Uc 

(w/m 2 .°C) 

1 - 8.67 

18.21 

30.37 

29 .16 

35.36 

33.72 

41 .41 

39.18 

41.19 

38.99 


Air with P- 100 psi (689 kpa) 
TABLE 3.2 -Continued 


m 

( Kg/Hr) 


0.6950 

0.9829 

1.2038 

1.3900 


(°C) 

52.40 

45.38 

44.62 

42.96 


(°C) 


(°C) 


I 


58.30 78.10 
42.20 46.71 
39.35 44.80 
36.45 39.78 



1.554 


40.00 


28.10 


34.50 


























TABLE 3.3 



1.5540 


41.34 


31.80 


36.28 




































TABLE 3.4 


m 

(Kg/Hr) 


0.6950 3070 


0.9829 


1 .2038 


1 . 3900 


1 .5540 


4341 


5317 


6139 


6864 


20.534 


38.08 


45.611 


51.28 


54.33 


Nu 

3. 

.091 

5. 

.731 

6, 

.865 

7. 

.718 

8 

.178 


18.99 


33.11 


36.67 


42 .67 


44.76 


Test #1, Air with P* 300 psi (2068 kpa) 


TABLE 3.4 -Continued 


UC 

( w/m 2 

•°C) 

18 . 

63 

32 . 

02 

37 . 

18 

40. 

86 


m 

Cg/Hr) 

Tc 

(°C) 

ATc 

(°c) 

ATm 

(°C) 

0.6950 

52 . 20 

55.82 

73.40 

0.9829 

49.40 

51.72 

51.87 

1.2038 

45.79 

41.33 

42.38 

1.3900 

44.20 

38.00 

40.02 

1.5540 

42.11 

34.00 

37.78 
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TABLE 3.5 




m 

(Kg/Hr) 

Re 

h 

( w/m 2 • °C ) 

Nu 

U 

(w/m 2 -°C) 

Uc 

( w/m 2 - °C ) 

0.6950 

3070 

24.596 

3.702 

22.42 

21 .98 

0.9829 

4341 

40.52 

6.098 

34.94 

33.88 

1 . 2038 

5317 

46.68 

7.027 

39.43 

38.08 

1.3900 

6139 

53.25 

8.016 

44.02 

42 . 35 

1 . 5400 

6864 

56.85 

8.557 

46.45 

44 . 59 


Test #1, Air with P- 400 psi (2758 kpa) 


TABLE 3.5 -Continued 


m 

(Kg/Hr) 

Tc 

(°C) 

ATc 

(°C) 

ATm 

(°C) 

0.6950 

48.40 

51.82 

56.89 

0.9829 

46.57 

51.14 

48.20 

1.2038 

42.04 

42.12 

42.20 

1.3900 

39.24 

33.25 

33.72 


1.5540 


38.14 


32.60 


34.62 








































TABLE 3.6 


m 

(Kg/Hr) 


0.6950 3070 


0.9829 4341 


1.2038 5317 


1.3900 6139 


1.5540 I 686 


( w/m^ • °C ) 


26.268 


42.03 


49.255 


55.103 


57.654 


3.954 


6.327 


7.414 


8.29 


8 


(w/m^*°C) (w/m^«°C) 


23.81 


36.06 


41 .25 


45.98 


46.98 


Test #1, Air with P» 500 psi (3448 kpa) 


TABLE 3.6 -Continued 


23.37 


35.07 


43 . 73 


45.33 


m 

(Kg/Hr) 

Tc 

(°c) 

4 Tc 
(°C) 

ATm 
( °C ) 

0.6950 

49.00 

50.12 

51.52 

0.9829 

47.10 

48.30 

43.88 

1.2038 

43.00 

42.40 

40.26 

1.3900 

40.81 

39.61 

38.82 

1.5540 

38.47 

34.95 

36.60 
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f# 


appendix c 

SAMPLE CALCULATIONS 
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Id 


I . Heat Transfer Analysis (Oil or Water Case) 
Data Available (from page 144 ) 

Re - 22 
P ■ 0 psi 
D ■ 0.00635 m 
q"- 2849.40 v/m’ 

K - 13.93 x 10~ 2 w/m .°C 
21 °C 
181.10 °C 


in 


T c - 85.60 °C 

Ax * 0.00787 m 

K . . „ » 2 v/m. °C 

graphite 


1- Heat Transfer Coefficient 


h - JL 


2849.40 


- 29.84 v/m .°C 


T w - T c (181.10 - 85.60) 


2- Nusselt Number 


Nu - M - (29. 84)(0. 00635) . i. 360 
K 13 x 10 -2 


3- Thermal Contact Resistance 




r c - 0.00357 e 


(-0.002356 P) 


r - 0.00357 °C.m’/v 
c 


/ 


4- 


Overall Heat Transfer Coefficient 


H + "T + ^ [ 2^T84 + °-003935 + 0.00357] 

u c - 24.36 v/m\°C 

f H * + 0 -003937] 

u - 26.70 w/m 1 . °C 


♦Similar procedure is repeated for water except for the use 
of the properties of water. 

11 • Heat Tra nsfer Analysis (Air Case) 

Data Available (from page 157 ) 

P ■ 100 psi 
Re - 4341 
m - 0.9829 Kg/hr 

AT m " T w - T c " 46 ‘ 71 ° C 

AT c * T out - T in » *2.20 °C 

A » 0.007187 m 3 ; D - 0.00558 m 

C p - 1.004 KJ/Kg. °C 

K a - 0.0371 w/m . °C 

Kg ■ 2 w/m.°C 

Ax - 0.00787 

Heat Transfer Coefficient 

h « (0.9829)(1. 004)(42. 20) [0 2781 
(0.007187)(46.71) 
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h ■ 34.50 v/m * . °C 

2- Nusselt Number 

NU . H . (0. 00558)04. 50) _ «. 1Q 
K 0.0371 

3- Thermal Contact Resistance 
r c - 0.0015812 e^ -0 * 00141 3 P ) 

r c - 0.0015812 °C.m’/w 

4- Overall Heat Transfer Coefficient 

U - l - 1 

c (i + A* . r ) ( 1 + 0.003935 + 

h K 34.5 

U c - 29.036 w/m J . °C 

U - 1 » 1 

(I + Ax) ( 1 + 0.003935) 

h K 34.5 

U ■ 30.37 w/m’^C 


0.0015812) 



